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EDITORIAL 


The need for candid and publicly available research in the natural sciences is as great as it has ever been. As 
I write this foreword, on my desk lies a report released a week ago by the Australian Academy of Science 
drawing attention to the low level of scientific literacy in our society — surely a warning signal given that 
numerous science- and environment-related issues are demanding thoughtful solutions. This week Professor 
Brian Schmidt of Mount Stromlo Observatory observed that research has a 10-to-20 year horizon and warned 
against “disjointed short-term initiatives” that destroy stability and confidence in the sector. 


Scientific research is a critical part of the infrastructure of a modern society, as much as the more readily 
recognised “hard” physical infrastructure such as roads, pipelines and ports. It forms the foundation of 
technological progress. It flashes amber lights about looming problems for the managers of land and sea. It 
supplies evidence vital for prudent policy formulation. As a public good, scientific research deserves to be 
funded generously from the public purse and scientists in the universities and public services warrant secure 
tenure and institutional support. 


So the papers in this latest Proceedings give cause for celebration. We have several interesting and content-rich 
papers which add substantially to the body of learned knowledge and also happen to have practical relevance 
for land managers and practitioners. The Society is proud to provide a vehicle for presenting these findings to 
Queensland and international audiences, continuing the tradition of its Proceedings for more than a century. 


The incoming Council’s primary new project in the forthcoming year is to develop a website which will serve 
as a repository for a range of scientific data, information and interpretative materials that can support scientists 
and decision-makers in the public and private sectors. A major objective of this initiative is to enhance general 
scientific literacy by improving public access to scientific evidence. 


Of course, another unwavering priority of the Society will be to publish Proceedings 119, for which I now 
invite contributions — both learned articles for peer review and opinion pieces on subjects relevant to the natural 
sciences in Queensland. 


Turning to procedural matters, the 2013 annual general meeting of the Society was hosted by the Queensland 
University of Technology at its Gardens Point campus. I thank Professor David Gust, Head of School — Earth, 
Environment and Biological Sciences, who discussed the role of societies like the Royal and their importance in 
research and in engaging a new generation of scientists. Prof. Gust then invited members and guests to tour the 
University’s impressive new Science and Engineering Centre. 


I pay aspecial tribute here to Craig Walton, outgoing President who for six years has discharged the responsibilities 
with good sense and good humour. Craig will continue on Council as Immediate Past President. Craig has 
ensured that the Society has maintained its status as the pre-eminent non-government scientific organisation in 
this State and has retained the respect of scientists. 


I also pay tribute to this issue’s Editor, Tanya Scharaschkin and the former editor, Raghu Sathyamurthy, who 
have guided the process through to production: science has been well served by their service. 


Geoff Edwards B.Sc.(Hons.); M.Pub.Ad.; PhD 
President 

Royal Society of Queensland 

August 2013 


Editors: S. Raghu! and T. Scharaschkin’? 
'Department of Entomology, Rice Research & Extension Center, University of Arkansas, Stuttgart, AR72160, USA 
?EEBS, Science and Engineering Faculty, Queensland University of Technology, Brisbane, OLD4001 


REGIONAL PANMIXIA IN THE MULLET MUGIL CEPHALUS 
ALONG THE COAST OF EASTERN QUEENSLAND; 


REVEALED USING SIX HIGHLY POLYMORPHIC MICROSATELLITE LOCI 


HUEY, J.A., ESPINOZA, T. & HUGHES, J.M. 


Catadromy describes an array of migration behaviours which are expected to generate different pat- 
terns of genetic structure; from total panmixia (genetic homogenisation) to moderate structure with iso- 
lation by distance. In this study, the genetic structure of the circumtropical sea mullet (Mugil cephalus) 
along approximately 550km of the Queensland coastline was investigated. Using six highly variable 
microsatellite loci, no significant genetic structure was detected (F’,,.< zero), despite high genetic diver- 
sity providing sufficient power to detect structure. These patterns of genetic structure were concordant 
with other studies of MZ cephalus at similar scales in other parts of the world. Overall, genetic structure 
for M. cephalus is similar to that of many marine species, which also show panmixia along this sec- 
tion of the Queensland coast. Conservation of this species should focus on ensuring freshwater-marine 
migration is not diminished. 


Joel A. Huey, Australian Rivers Institute, Griffith School of Environment, Griffith University. 
Nathan, Queensland, Australia, 4111, email: j.huey@griffith.edu.au; Thomas Espinoza, Queensland 
Department of Environment and Natural Resource Management. Bundaberg, Queensland, Australia, 
4670; Jane M. Hughes, Australian Rivers Institute, Griffith School of Environment, Griffith University. 


Nathan, Queensland, Australia. 4111 


INTRODUCTION 

Diadromous species require marine and freshwater 
environments to complete key components of their 
life history. For the conservation and management 
of such species, the quality of marine and freshwater 
environments is important, as is the maintenance of 
connectivity between these environments (McDowall, 
2008). When freshwater connectivity is reduced, 
saltwater connectivity among populations along 
a coastline becomes an important factor for long 
term population viability. However, despite the 
global distribution of diadromous species, and their 
prominence in freshwater and near shore ecosystems, 
little is known about the hydrological and genetic 
connectivity of populations (McDowall, 1999). 


Broadly, diadromy consists of three types; anadromy, 
amphidromy and catadromy. Catadromous species 
are the rarest of the three diadromous categories, with 
breeding occurring in marine environments followed 
by migration into rivers and lakes for adult growth 
(McDowall, 1992). The most extreme catadromous 
species, the freshwater eel (Anguilla spp), migrates 
to a single location in the ocean to spawn before 
recruits return to freshwater (Schmidt, 1923). The 
result is panmixia (genetic homogenisation) over the 
entire distribution of the species, across thousands 
of kilometres (Avise et al., 1986; Dannewitz et al., 


2005). However, not all catadromous species exhibit 
such extreme migratory behaviour, with many species 
Spawning in near-shore and estuarine environments. 
Modest levels of genetic structure ensue, with partially 
isolated populations exchanging some migrants via 
the marine environment. This leads to restricted gene 
flow and isolation by distance; a correlation between 
genetic and geographic distance (Chenoweth & 
Hughes, 1997; Jerry & Baverstock, 1998). Interpreting 
genetic pattern and process in catadromous species 
is evidently complex, yet essential to the effective 
management and conservation of these species. 

Migratory behaviour of the widely distributed sea 
mullet (Mugil cephalus L. 1758) was noted by the 
first natural historian, Aristotle who compared it to 
the freshwater eel (Anguilla anguilla L. 1758) in his 
Historia Animalium. A tropical and temperate species, 
M. cephalus is restricted to coastlines between 42 
degrees north and south latitudes (Thomson, 1963). 
Mugil cephalus \arvae hatch at sea and use tidal 
currents to enter estuaries, both passively and actively, 
at 28-42 days old (Pusey et al., 2004). Movement 
towards more turbid freshwater habitats is believed 
to be a result of osmoregulatory, thermal and trophic 
factors (Oren, 1981; Thomson, 1955). The majority 
of juvenile mullet will then remain within one river 
system until maturity. However, migrations are 
common, including ‘hard gut’ runs (also know as 
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‘wash out’ runs; believed to be a response to flood- 
induced removal of detrital food sources from riverine 
habitat) which may lead to colonisation of new river 
systems (Hadwen & Arthington, 2007; Pusey et al., 
2004; Smith & Deguara, 2002). When adults migrate 
back to the marine environment to spawn, movement 
is almost always in a northerly direction, with 
resulting spawning aggregations close to river mouths 
(Kesteven, 1953; Virgona et al., 1998). With all of 
these opportunities for gene flow among populations, 
it is likely that genetic structure among populations 
will be low to non-existent over extensive geographic 
areas joined by continuous coastlines. 


To date, population genetic studies of MZ. cephalus 
have used allozymes, mtDNA and AFLPs to detect 
genetic structure. A recently available microsatellite 
library for M cephalus (Miggiano et al., 2005) 
provides the first opportunity to explore fine-scale 
patterns of genetic structure for this species with 
microsatellites to more precisely test the hypothesis 
that populations are panmictic along coastlines. 


MATERIALS AND METHODS 
SAMPLE COLLECTION 
A mixture of adults and juveniles were collected from two 
sites in Baffle Ck (14 from the estuarine Ferry’s Crossing; 
and 6 from the upstream limit; combined 1 = 20), one 
reach in the lower Boyne (fresh, 7 = 21), one site Calliope 


Queensland 


O Study Site 


New South Wales 
Approx 200km 





FIG. 1. Map of study area showing sites. Site latitude and 
longitudes are provided in Table 1. Grey area on the map of 
Australia (inset) represents the distribution of 4 cephalus 
in Australia. 


(first fresh waterhole above the estuary 7 = 8), one site in 
the Albert River (freshwater 7 = 18) and one site in the 
Logan River (freshwater n = 20) (Fig. 1). Samples were 
taken from riverine sites as M. cephalus was aggregated 
into discrete locations at this life-history phase. 


Sampling was undertaken using a boat electrofisher 
(2.5GPP, Smith-Root, Inc. Vancouver, WA, USA) 
between May 2009 and May 2010. Fish were held 
within an on-board aerated tank (200L) fitted with 
a recirculation system (Flow-Rite Controls, Ltd. 
Australia) before being anaesthetised using AQUI-S 
(540g/L Isoeugenol) solution (Sml/1000L). Once 
sedated a 5xSmm sample was taken from the anal fin 
using sterilised scissors and placed in a labelled vial 
containing absolute-ethanol (C,H.OH = 46.1). 


LABORATORY METHODS 

Total genomic DNA was extracted using a salt 
extraction technique, adapted from Aljanabi and 
Martinez (1997). After precipitation and washing, 
DNA was resuspended in 50uL ddH,O until required 
for further analysis. 


Eight microsatellite loci were taken from Miggiano 
et al. (2005) and pre-screened to identify which loci 
amplified in the samples. Loci 15CM, 16DM, 17FM, 
2DM, 2FH and 6DM were considered appropriate 
for this study, displaying high polymorphism and 
easily scored alleles. Loci ISAM and 1EH failed 
to amplify successfully. Screening of microsatellite 
loci was performed on an Applied Biosystems 
3130 Genetic Analyser; which requires amplified 
fragments to be fluorescently labelled. To achieve 
this, forward primers for each locus had a 20- 
mer oligonucleotide chain added, which, with the 
addition of a fluorescently labelled tag in the PCR, 
enabled the amplified fragment to be fluorescently 
labelled. A different 20-mer oligonucleotide chain 
was used for each dye (6-FAM, VIC, NED, PET) as 
outlined in Real et al. (2009). Multiplex PCR’s were 
not used to amplify loci. 


10ul PCR = reactions included; O.5uL template 
DNA, 0.5mM reverse primer, 0.125mM_ forward 
‘tailed’ primer, 0.4mM _ fluorescently labelled tag, 
0.2mM dNTP’s, 1.5mM MgCl, 1x Buffer, 0.2ul 
taq polymerase (red tag, Astral Scientific). Cycling 
conditions consisted of an initial hold at 94°C for 
5 minutes, followed by 40 cycles of 94°C for 30 
seconds, 50°C for 30 seconds and 72°C for 45 
seconds. There was a final extension for 7 minutes, 
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TABLE 1. Study site details and samples sizes 
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followed by an infinite hold at 4°C. PCR product was 
kept at 4°C until analysed. Loci were screened on an 
Applied Biosystems 3130 Genetic Analyser, with all 
loci run simultaneously. Genotypes were analysed and 
exported using GENEMAPPER (Applied Biosystems, 
version 4.0). All individuals were screened at all loci. 


STATISTICAL ANALYSIS 

The microsatellite dataset was tested for null 
alleles, poor scoring and large allele dropout using 
MICROCHECKER (van Oosterhout et al., 2004). 
Genetic diversity in the form of observed and expected 
heterozygosity was calculated using ARLEQUIN 
(version 3.5.1.2, Excoffier & Lischer, 2010). Allelic 
Richness was calculated in FSTAT (version 2.9.3.2, 
Goudet, 1995) with the estimate scaled by the smallest 
site’s sample size (Calliope, 7 = 8). Deviations from 
Hardy-Weinberg Equilibrium (HWE) were tested 
in ARLEQUIN using exact tests. The inbreeding 
coefficient F,, was calculated in GENEPOP (Raymond 
& Rousset, 1995; Rousset, 2008). 


Genetic structure was estimated using Weir and 
Cockerham’s 0 (Weir & Cockerham, 1984), an F,, 
analogue performed in ARLEQUIN. For the rest of 
the manuscript 0 will be referred to as F.,. Overall 
and pairwise F, values were calculated. Statistical 
significance of the F’,,. estimate was calculated using 
permutation tests. Furthermore, to identify significant 
structure among sites, Fisher’s exact tests were 
calculated using ARLEQUIN. 


As very low genetic structure is predicted for M 
cephalus across the study area it is important to 
delineate between insufficient statistical power 
revealing insignificant genetic structure and genuine 
panmixia. To do this a post-hoc power test was 
performed in POWSIM (Ryman & Palm, 2006). 
POWSIM generates many datasets with specific 


F,, values and then identifies what proportion of 
times, with a certain configuration of data (number 
of alleles, allelic frequencies, number of loci) 
significant structure will be detected at a specified 
confidence level (a = 0.05). The F,, values are 
derived by altering key parameters in a coalescent 
model; N., the effective population size, and t, the 
time of population divergence. 


For this analysis, different F’,, values were chosen 
to identify the point at which the data became 
non-informative (low power). Allele frequency 
distributions, numbers of alleles, sample sizes 
and numbers of loci for each simulation matched 
those found in this study allowing an estimation 
of power specific to our data configuration. To 
generate simulated data sets with these F,, values, 
the equation outlined in Ryman and Palm (2006) was 
used (Equation 1). 


F,, = 1-(1-(1/2N,)) 


EQUATION 1 


Multiple configurations of t (time since divergence) 
and N_. (effective population size) can generate the 
same Ff’... Running the analysis with a smaller N, 
makes the analysis quicker, but may lead to the loss 
of low frequency alleles (Ryman & Palm, 2006). 
For this analysis N, was set at 10,000, ensuring the 
retention of low frequency alleles. It is important 
to note that the N, used for this analysis does not 
need to match expectations for wild populations. 
As N, was fixed (10,000) for each power test 
(different F,, values), only t was altered to generate 
the appropriate FP. For each test, 1000 sample 
data sets were generated and analysed using F,,. 
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x? permutation tests and Fisher exact tests. These 
are both calculated within POWSIM. The relative 
frequency of significant tests was determined and 
used as a close approximation of the power for the 
data to detect structure at that F,, value. 


RESULTS 

GENETIC DIVERSITY 

All individuals amplified successfully at every locus. 
MICROCHECKER uncovered no evidence for null 
alleles, poor scoring or long allele drop out. As such, 
all loci were included in all analyses. Genetic diversity 
was high across all loci and populations (Table 2). 
Expected heterozygosities ranged between 0.87 and 
0.97, with allelic richness (scaled by the smallest site; 
Calliope, 7 = 8) ranging between 8.3 and 12.1 alleles. 
These levels of diversity were broadly consistent with 
those detected in the primer note (Miggiano et al., 
2005). Two populations deviated significantly from 
HWE ata locus. The site ‘Logan’ exhibited a deviation 
at locus 17FM, with this population exhibiting a 
deficiency of heterozygotes (Ff, = 0.15). Miggiano 
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et al. (2005) also detected a heterozygote deficiency 
with this locus in the Sardinian population, but not 
a population from Moreton Bay. The site ‘Baffle’ 
deviated from HWE at locus 2FH. 


GENETIC STRUCTURE 

F,,, permutation tests and Fisher’s exact tests revealed 
no significant genetic structure (P < 0.05). All Fy, 
values (overall and pairwise) were negative. While 
F,,, is generally accepted to range between 0 and 1 (0 
= genetic homogenisation), ARLEQUIN can return 
values less than zero in some circumstances. These 
values should be interpreted as zero (Meirmans, 2006). 


STATISTICAL POWER 

Using specific F,, values, POWSIM described an 
exponential relationship between F,,, and statistical 
power (Fig. 2). Both methods for statistically 
detecting structure, y? permutation tests and Fisher’s 
exact tests, revealed the same distribution. At very 
small F’... values (e.g., F’,,, between 0 and 0.0005) the 
proportion of significant tests closely mirrored that 





Fig. 2. Power analysis results. Line describes the relationship between F’, and the proportion of significant tests using the 7 
contingency test. The solid circles describe statistical power using the observed data to parameterise the analysis. The hollow 
circles describe statistical power with all sample sizes increased to 50 individuals. 
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Table 2. Diversity indices showing expected heterozygosity (H,,), Allelic Richness (rarefied to n=8, A,), the inbreeding 
coefficient (F x) and HWE p-values 
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HWE p-values 





expected by the Type I error rate of 0.05 (That is, 5% 
of tests should be significant in the absence of true 
differentiation with an alpha value of 0.05). For this 
data, population structure can be detected at sufficient 
power using six microsatellite loci when F,, values 
are above 0.008 (proportion of significant tests = 
0.934 and 0.927; for y° and exact tests respectively). 


As the sample sizes for this study were lower than 
what is generally collected for gene flow studies (e.g., 
typical studies of this nature may aim for >30 samples 
per study site), we also performed a power analysis 
with the sample sizes increased to 50 individuals per 
population. The allelic distributions were kept the 
same and the analysis was only conducted using the 
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y° test. As expected, greater power resulted (Fig. 2). 
More than doubling the sample sizes allowed reliable 
detection of genetic structure at /’,,, values as low as 
0.003 (proportion of significant tests = 0.973) 


DISCUSSION 

Mugil cephalus, despite breeding in semi-discrete 
locations along the coast, revealed no evidence 
of genetic structure across ~550 km of coastline. 
This was despite the use of six highly variable 
microsatellite loci. A post-hoc analysis of power 
revealed that with these loci and sample sizes, 
genetic structure could be reliably detected at F,,. 
values above 0.008 (proportion of significant tests = 
0.934 using x? tests). As such, it is clear that for this 
data, only very weak levels of genetic structure (< 
0.008) would go undetected (1.e., be non-significant). 
A further power analysis showed that detectable 
structure could be revealed with F,,, values as low as 
0.003, if the sample sizes were more than doubled (n 
= 50 per site). While it is inevitable that population 
genetic studies will always have some small levels 
of genetic structure that will go undetected due to 
insufficient power, it is important to consider the 
biological significance of these very low F,, values. 
For example, in the truly panmictic European eel 
(A. anguilla) microsatellite loci have revealed 
significant F,, values exceeding 0.008 (Dannewitz 
et al., 2005; Wirth & Bernatchez, 2001). These 
significant and small fF’, values (as low as 0.0006 
and still significantly different from zero due to high 
statistical power) were not the product of restricted 
gene flow (as proposed by Wirth & Bernatchez, 
2001), but were the product of temporal instability in 
gene frequencies (Dannewitz et al., 2005). 


Previous studies of M. cephalus have found 
conflicting patterns of genetic structure. Along the 
east-coast of China, studies have revealed some 
evidence for restricted gene flow with isolation by 
distance (Jamandre et al., 2009; Lui et al., 2009; 
Lui et al., 2009). These results differ from the 
findings of this study, and from other population 
genetic studies of M. cephalus which have found 
no genetic structure (Rocha-Olivares et al., 2000). 
However, considering that M cephalus may 
represent a morphologically conserved complex of 
cryptic species (Crosetti et al., 1994; Jamandre et 
al., 2009; Livi et al., 2011), differences in spatial 
genetic structure for different regions may be due 
to different species, which might have different 
dispersal capacities and behaviours. 


Other studies involving catadromous species along 
the Australian coast have revealed significant 
genetic structure at comparable geographic scales. 
Furthermore, these studies have generally revealed 
evidence for isolation by distance. For example, along 
the same coastline, studies of the Australian bass 
(Macquaria novemaculeata Steindachner 1866) have 
identified restricted gene flow using allozymes and 
mtDNA sequence data (Chenoweth & Hughes, 1997). 


When considering genetic structure in aquatic 
organisms, it is generally accepted that freshwater 
species tend to be highly structured compared to 
marine species (de Woody & Avise, 2000; Ward et 
al., 1994). Diadromous species exhibit both patterns, 
ranging from the highly structured (e.g., anadromous 
salmon) to the entirely panmictic (e.g., catadromous 
anguillid eels). The patterns of genetic structure 
found for M. cephalus along the eastern coast more 
closely match those of marine species, with very 
low levels of structure detected over vast geographic 
areas. Marine species generally exhibit no, or weak, 
genetic structure along this section of the east 
Australian coast. This includes an array of species 
and dispersal strategies; including, fish (Sumpton et 
al., 2008), crustaceans (Gopurenko & Hughes, 2002; 
Haig et al., 2010), and bivalves (Murray-Jones & 
Ayre, 1997). This concordance between species with 
such an array of life history strategies suggests that 
currents in this area are sufficient to facilitate gene 
flow up and down this section of coast. However, for 
the fish (snapper, Pagrus auratus Bloch & Schneider 
1801) and crustaceans (Phycomenes zostericola 
Bruce 2008 and Scylla serrata Forskal 1775), broader 
geographic scales revealed considerable genetic 
structuring (Gopurenko & Hughes, 2002; Haig et al., 
2010; Sumpton et al., 2008). Further research on M. 
cephalus stocks in central and northern Queensland, 
together with NSW are warranted to determine if the 
lack of genetic structure found in this study exists at 
larger spatial scales. 


The results of this study have important implications 
for informing our understanding of the biology of M. 
cephalus. Firstly, it suggests that both ‘hard gut’ and 
spawning migrations undertaken by juvenile and adult 
sea mullet are sufficient to not only homogenise gene 
frequencies at significant spatial scales (~550km), 
but also subsidise genetic vigour into adjacent river 
systems with reduced riverine-marine connectivity. 
The northerly direction of these volitional migrations 
(up to 700km; Pusey et al., 2004) may be equally 
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reciprocated by the southerly impetus of inshore 
currents along the east Australian coast. Importantly, 
while connectivity up and down the coast may be 
important for population viability in M. cephalus, 
connectivity between marine and freshwater habitats 
may be the most important process for this migratory 
Species. For example, the installation of water 
infrastructure in riverine and estuarine habitats is 
ongoing in eastern Australia and poses an immediate 
threat to diadromous species in Australia (Crook et al., 
2006; McDowall, 1999; Schmidt et al., 2011; Stuart & 
Berghuis, 2002). Indeed, evidence already exists for 
the deterioration of freshwater stocks of MZ cephalus 
in Queensland due to artificial in-stream barriers 
(Heidenreich & Broadfoot, 2001; Stuart & Berghuis, 
1997; 1999). Therefore, while locally extinct rivers 
have the potential to be recolonised, the persistent 
impact of in-stream barriers may continually inhibit 
the ability for migrating individuals to achieve this. 
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DO CANE TOADS (RHINELLA MARINA) 
IMPACT DESERT SPRING ECOSYSTEMS? 


CLIFFORD, S.E., STEWARD, A.L., NEGUS, P.M., BLESSING, J.J. 
& MARSHALL, J.C. 


Since their introduction in 1935, cane toads (Rhinella marina (Linnaeus, 1758)) have established and 
spread throughout north and north-eastern Australia. Cane toad impacts to terrestrial ecosystems are 
well documented, but impacts to aquatic ecosystems are less well known. We investigated the diet of 
cane toads collected from warm Great Artesian Basin-fed springs on Edgbaston Reserve in Central 
Queensland, Australia. A higher proportion of aquatic invertebrates to terrestrial invertebrates were 
found amongst their alimentary canal contents. Aquatic taxa consumed included molluscs (Gastropo- 
da), insects (Coleoptera) and crustaceans (Amphipoda). Given this diet, the presence of cane toads at 
Edgbaston Springs, and the high endemicity of the aquatic biota of these springs, we conclude that R. 
marina present a threat to the conservation of desert spring ecosystems. 


Sara E. Clifford (sara.clifford@qld.gov.au), Alisha L. Steward, Peter M. Negus, Joanna J. Blessing, 
Jonathan C. Marshall; Queensland Department of Science, Information Technology, Innovation and 
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the Arts, GPO Box 5078, Brisbane, OLD, 4001. 


INTRODUCTION 

The cane toad Rhinella marina (formerly Bufo 
marinus (Linnaeus, 1758)) is an amphibian native to 
Central and South America belonging to the family 
Bufonidae. This species was introduced to coastal 
Queensland in 1935 as an ultimately unsuccessful 
biological control agent for sugar cane pests (Freeland, 
1984; Lever, 2001). Since then, it has spread widely 
throughout north and north-eastern Australia (Sutherst 
et al., 1995; Urban et al., 2007) causing significant 
negative impacts to Australian ecosystems (Phillips et 
al., 2003; Shine, 2010). 


Cane toads are opportunistic generalist feeders (Zug 
& Zug, 1979; Reed et al., 2007; Heise-Pavolv & 
Longway, 2011) and are able to withstand a wide 
range of climatic conditions (Lever, 2001; Urban et 
al., 2007). Their resilient nature in combination with 
the high vagility of adults has enabled the species to 
expand its range to now occupy over 1.2 million km? 
of Australia. They are predominantly found in tropical 
and subtropical areas including much of Queensland. 
There is potential for this range to further expand to 
over 2 million km’ across all mainland states (Urban 
et al., 2007). 


Like many introduced species, cane toad populations in 
Australia are exposed to few of the predators, parasites 
and pathogens present in their native range (Speare, 
1990), reducing the impacts of predation and disease 
on population size and life expectancy. Furthermore, 


all life stages of the cane toad (from egg to adult) 
contain toxins such as bufotoxins and bufogenins, 
which deter predators (Alford et al., 1995). Having 
not historically encountered these toxins, Australian 
fauna lack behavioural or physiological mechanisms 
to alleviate either exposure or the toxic responses to 
exposure (Lever, 2001). 


While the ingestion of cane toad toxins threatens a 
variety of Australian fauna (Doody et al., 2009; Phillips 
et al., 2003; Shine, 2010), the direct consumption of 
native fauna by the cane toads is also a threat. Cane 
toads are indiscriminate feeders able to form large 
populations, which can consume large quantities of 
invertebrates (Zug & Zug, 1979; Freeland et al., 1986; 
Shine, 2010; Heise-Pavlov & Longway, 2011). 


The adult cane toad diet has been studied extensively, 
with alimentary canal contents dominated by 
terrestrial invertebrate prey (e.g. beetles, termites, 
ants) (Freeland, 1984; Striissmann et al., 1984; 
Freeland et al., 1986; Reed et al., 2007; see also Shine, 
2010). Current evidence is lacking regarding adult 
cane toads specifically targeting aquatic invertebrates, 
with no previous research into the diet of cane toads 
inhabiting springs. There are however documented 
cases of aquatic macroinvertebrates, such as beetles 
from the families Hydrophilidae and Dytiscidae, being 
consumed by cane toads (Hinckley, 1963), indicating 
this species is a potential threat to aquatic ecosystems. 
Aquatic predators are further threatened by cane toads 
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FIG. 1. Cane toads and macroinvertebrates were sampled 
from within a Great Artesian Basin spring wetland located 
on Edgbaston Reserve, Queensland. 


as the consumption of cane toad eggs and tadpoles 
can be fatal due to their toxicity (Crossland & Alford, 
1998; Somaweera & Shine, 2012). 


R. marina was recently identified as a potential 
threat to the unique Great Artesian Basin (GAB) 
spring wetland communities in Queensland 
(Fensham et al., 2010), including the Edgbaston 
spring complex. The spring wetlands within 
Edgbaston Reserve represent a permanent source of 
water in a semi-arid region, and their long isolation 
has resulted in a diverse and endemic GAB spring 
community of fish, aquatic plants and aquatic 
macroinvertebrates (Ponder & Clarke, 1990; 
Wager & Unmack, 2000; Fensham et al., 2010). All 
GAB spring communities are listed as endangered 
under the Australian Environment Protection and 
Biodiversity Conservation (EPBC) Act 1999, 
with several endemic species also individually 
listed as endangered or vulnerable (of particular 
note are red-finned blue-eye (Scaturiginichthys 
vermeilipinnis) and the Elizabeth Springs goby 
(Chlamydogobius micropterus)). 


As available sources of water in often arid or semi-arid 
regions, GAB springs inherently attract a range of fauna 
including introduced species such as feral pigs and cane 
toads (Fensham et al., 2010; Fensham et al., 2011). 
While the permanency of surface water in springs 
provides opportunity for cane toads to persist in this dry 
landscape, the underlying processes and mechanisms of 
their potential impacts are yet to be quantified. 


This paper presents an initial investigation of the 
adult cane toad diet within a GAB spring wetland 
at Edgbaston Reserve in Central Queensland. It was 
hypothesised that given the opportunistic feeding 
habits of adult cane toads and the limited surface 
water in these regions, adult toads will be consuming 
aquatic invertebrates from GAB springs. If so, 
this poses a threat to the endangered GAB spring 
community within Edgbaston Reserve, in particular 
the rare and endemic aquatic invertebrate taxa. 


MATERIALS AND METHODS 

STUDY LOCATION 

Toads were collected from one GAB spring located 
within Edgbaston Reserve (Figure 1) approximately 80 
km northeast of Aramac, in Central Queensland, Australia 
(22°44’°E, 145°25’S) (Wager & Unmack, 2004). The 
Edgbaston spring complex is on the eastern side of the 
Reserve at the base of the Desert Uplands escarpment, 
on the north-eastern side of the GAB (Kerezsy, 2011). 
There are up to 180 springs in the complex (Fensham & 
Fairfax, 2009) with varying surface extent. These contain 
the highest number of endemic macroinvertebrates of all 
the spring complexes in Australia (Ponder et al., 2010) and 
are home to two endemic fish species — the endangered 
red-finned blue-eye (S. vermeilipinnis) and the vulnerable 
Edgbaston goby (C. squamigenus). 


Cane toads and macroinvertebrates were collected from 
spring NW30, which is situated in close proximity to other 
springs in the northern section of the Edgbaston spring 
complex. NW30is one of the larger springs in the complex, 
with a surface extent of 2,723 m?at the time of sampling 
(mean extent within complex = 1,155 m’) (Blessing et al., 
2012). Harsh weather conditions characterise the semi- 
arid climate in this area. Mean annual rainfall is 692 mm 
while mean annual evaporation is over 1,997 mm (EHP, 
2009). Mean minimum and maximum temperatures in 
the area (Barcaldine) for the month of sampling (July) are 
7.9 °C and 22.6 °C respectively (Bureau of Meteorology, 
2012). Spring NW30 is fed by warm groundwater and has 
a daytime water temperature range of 20-33 °C throughout 
the year (Fairfax et al., 2007). 
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FIG. 2. Photographs of a Great Artesian Basin spring wetland within Edgbaston Reserve showing (A, B) a warm spring vent; 
(C) spring vent encompassing the area sampled for toads and aquatic macroinvertebrates; (D) cooler spring outflow area — no 
toads were captured in this area. 


SAMPLE COLLECTION 

Cane Toads 

Cane toad specimens were collected during a two- 
person 15 minute spotlighting survey in and around 
the entire extent of spring NW30 (Figure 2), during 
the late evening in early July, 2011. All specimens 
were euthanised and frozen prior to transportation to 
the laboratory. 


Aquatic Macroinvertebrates 

Aquatic macroinvertebrates were collected from 
spring NW30 using a 250 um mesh dip net. Five 
metres of spring habitat were sampled using a 
combination of short lateral sweeps (approximately 
30 cm each) and vertical lifts. Macroinvertebrates 
were live-picked in the field and preserved in 
100 % ethanol for transportation. Specimens 


were typically identified to the taxonomic level 
of family, with the exception of Chironomidae 
(identified to subfamily), Acarina, Hirudinea, and 
Oligochaeta (identified to subclass) and Ostracoda 
(identified to class). 


LABORATORY ANALYSIS OF CANE TOADS 
Cane toads were defrosted prior to laboratory analysis, 
blotted dry with paper, weighed to the nearest 0.1 
g and snout-urostyle length (SUL) (mm) recorded 
prior to dissection. The length of the alimentary canal 
(mouth to anus) was removed, measured and placed 
onto a Petri dish lined with 1 mm graph paper. Visual 
assessments of both the stomach and intestinal tract 
were made to determine the fullness of each using the 
following categories: 1) Empty; 2) Little (<25 %); 3) 
Some (25-75 %); and 4) Full (©75 %). 
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An indirect volumetric method was used to assess the 
relative contribution of food items in the alimentary 
canal as per Hyslop (1980). Alimentary canal contents 
were placed onto a Petri dish sitting on top of | mm 
graph paper and compressed to a constant depth 
of approximately 1 mm thickness. Each food item 
was scored according to the number of graph paper 
squares covered, and expressed as a percentage of the 
total number of squares covered. The volumes of large 
items that could not be compressed were estimated. 


Food items were categorised into: unidentified 
material, detritus and sand, aquatic invertebrates and 
terrestrial invertebrates. Where possible, individual 
Specimens were further identified to class, order or 
family level and the number of each noted. Coleoptera 
were identified to family where possible; however, in 
the case of detached elytra, these were classified as 
‘aquatic’ if they were visually identical to the elytra of 
beetles in the corresponding aquatic macroinvertebrate 
sample or to other identified aquatic specimens within 
the toad alimentary canal contents. Two detached 
elytra were counted as representing a whole specimen. 


RESULTS 

CANE TOADS 

Thirteen cane toads (12 male and 1 female) were 
collected from spring NW30 during the survey, all of 
which were found within three metres of the spring 
vent despite the total surveyed area being larger. The 
mean defrosted specimen weight was 45 g (range 
24.7 g — 134.3 g); mean SUL was 76 mm (range 64 
mm — 111 mm): and the mean alimentary canal length 
282 mm (range 198 mm — 414 mm). Three toads had 
empty stomachs, while the remaining 10 contained 
‘little’ stomach matter (i.e. <25 % full). Intestinal tract 
fullness consisted of 4 toads with ‘little’ content and 9 
with ‘some’ content (1.e. 25-75 % full). 


Of material identified from cane toad alimentary canal 
contents, 65 % of the volume fell into the category of 
‘detritus and sand’ (Figure 3). Less than 2 % of the material 
was categorised as ‘undetermined digested material’ that 
could not be identified, and this was recorded in only three 
of the toads. The remaining 34 % of the alimentary canal 
contents was identified as invertebrates represented by 
fifteen aquatic and four terrestrial taxa. Of this volume, 
aquatic invertibrate taxa made up 29.5 %, whereas 
terrestrial invertebrates made up only 4 %. Aquatic 
Coleoptera represented 15.8 % of the volume, aquatic 
Gastropoda 11.4 %, and the remaining 2.3 % of aquatic 
invertebrates consumed was comprised of Acarina, 
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FIG. 3. Mean percentage of food items found in cane toad 
alimentary canal contents (v = 13). Standard error bars (+ 1 
standard error) are shown. 


Amphipoda, Diptera, Epiproctophora, Hemiptera, 
Hirudinea, and Oligochaeta (Table 1). On average, toad 
alimentary canal contents contained approximately 40 
individual aquatic organisms (dominated by Coleoptera 
and Gastropoda) and 2 individual terrestrial organisms 
(dominated by the family Formicidae i.e. ants) (Table 1). 


AQUATIC MACROINVERTEBRATES 

Fourteen taxa were represented in the aquatic 
macroinvertebrate sample collected from 
Spring NW30, (Table 1). Specimens from the 
families Dugesiidae, Leptoceridae, Culicidae 
and from the class Ostracoda were present in the 
macroinvertebrate sample but not found within the 
alimentary canal contents. Conversely, terrestrial 
invertebrates and several aquatic macroinvertebrates 
were found exclusively in the cane toad alimentary 
canal contents. 


DISCUSSION 

Previous studies have identified _ terrestrial 
invertebrates as the dominant prey items of cane toads, 
though as opportunistic feeders they have the ability 
to significantly impact other taxa that would normally 
comprise only a small part of their diet (Shine, 2010). 
Alimentary canal contents from this study confirmed 
this as, with the exception of detritus and sand, 
aquatic invertebrates made the largest contribution to 
the diets of R. marina collected from a GAB spring 
within Edgbaston Reserve. This was in terms of both 
contributions to alimentary canal volume and the 
number of individual prey items consumed. 


Cane toads were found to consume a large proportion 
of the available aquatic taxa, with eight of the 
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Table 1. Aquatic and terrestrial invertebrate taxa recorded from cane toad alimentary canal contents and the aquatic invertebrate 
sample collected from the corresponding GAB spring. Mean abundances of invertebrate taxa consumed per cane toad are 
shown, with ranges (minimum and maximum) in brackets. 
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eleven orders of aquatic invertebrates found in the 
corresponding macroinvertebrate sample also found 
in alimentary canal contents. Previous studies have 
shown that toads will preferentially consume small- 
bodied terrestrial invertebrates (Strussmann et al., 
1984): however, they have also been shown to have 
opportunistic generalist feeding habits (Zug & Zug, 
1979; Strussmann et al., 1984; Reed et al., 2007; 
Heise-Pavolv & Longway, 2011) which these results 
support. The categorisation of over half the alimentary 
canal contents as ‘detritus and sand’ (mainly sand) 
demonstrates the ‘sloppy’ feeding style of toads (Zug 
& Zug, 1979), which results in them often ingesting 
large amounts of superfluous material such as sand 
and detritus when capturing their prey (Hinckley, 
1963; Zug & Zug, 1979). 


Aquatic beetles (Coleoptera) and snails (Gastropoda) 
accounted for the majority, by volume and number, 
of the aquatic taxa consumed. Of the toads examined, 
85 % had aquatic beetles in their alimentary canal 
contents. Although gastropods were not identified to 
species level for confirmation of endemicity (due in 
part to damage during digestion), it is likely endemic 
Species were consumed since all six snails within the 
family Hydrobiidae found at Edgbaston Springs are 
endemic, as well as one species of Bithyniidae and 
several species of Planorbidae (Ponder et al., 2010). 


The distribution of material in the alimentary canal 
showed digestion of prey had already progressed 
beyond the stomach to the intestine, meaning soft 
bodied taxa such as flatworms (the only planarian 
found at Edgbaston is the endemic Dugesia artesiana 
(Sluys et al., 2007)) could potentially have been 
consumed but already digested beyond identification 
(Heise-Pavlov & Longway, 2011; Zug & Zug, 1979). 
To minimise this, future studies could consider instant 
freezing using liquid nitrogen or the in-situ removal 
of the alimentary canal prior to transportation to the 
laboratory for processing. 


If the abundance of aquatic macroinvertebrates 
consumed by cane toads continues at this rate, 
there is potential for cane toads to alter the spring’s 
macroinvertebrate community. The significance of 
this finding is compounded given that the collection of 
toads for this study was undertaken during mid-winter 
when toad activity is suppressed by cold nightly air 
temperatures (Freeland, 1984) (below 0 °C at the time 
of collection). It is of interest to note that although 
the sampling area encompassed the entire spring 


extent, cane toads were only found and collected in 
close proximity to the spring vent. The temperature 
of the discharging groundwater was warm (24 °C) and 
remains above 20 °C throughout winter (Fairfax et al., 
2007). With increasing distance from the spring vent, 
water temperatures dropped significantly and no toads 
were detected more than three metres from the warm 
vents. During warmer ambient conditions or periods of 
higher rainfall, toads may change their diets (possibly 
to more terrestrial sources lessening the pressure on 
the springs) as their ability to feed away from the 
warm spring water increases. Further diet analyses 
of additional toads collected from in and around the 
springs during warmer months and in periods of 
high rainfall, along with concurrent terrestrial and 
aquatic invertebrate sampling, are required. This will 
reveal if cane toads continue to preferentially feed on 
aquatic macroinvertebrates throughout the warmer/ 
wetter months, exerting ongoing pressure on endemic 
species and communities. 


Modelling using biophysical and climatic data shows 
that much of Queensland and Northern Australia is 
currently suitable for cane toads, and will continue to 
be suitable under future climate scenarios (Kearney 
et al., 2008). Currently found within Eastern GAB 
springs in Queensland (Fensham et al., 2010), cane 
toads have the potential to continue expanding their 
range to other GAB spring communities and wetlands 
in regions characterised by a scarcity of surface water. 


This initial investigation supports the notion that cane 
toads can directly impact GAB spring communities 
via predation of aquatic invertebrates. It is possible 
that the local consequences of this could be significant, 
given the small spring size and the endemicity of the 
aquatic invertebrate fauna. The trophic cascade caused 
by the feeding habits of cane toads could also pose a 
threat to spring communities, as experimental studies 
have attributed changes to terrestrial/floodplain 
invertebrate community assemblages to cane toad 
predation (Greenlees et al., 2006; Shine, 2010). 
Large numbers of cane toads in areas of the Northern 
Territory have also been identified as the apparent 
cause of reductions to both the abundance and species 
diversity of insectivorous reptiles due to the depletion 
of their food supply (Catling et al., 1999). 


These results suggest that in addition to current 
practices to manage threats (see Kerezsy, 2011), 
management of the Edgbaston spring complex could 
consider incorporating measures to reduce toad 
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abundance or prevent them from accessing springs. 
That cane toads seem to congregate around the spring 
vents could also aid in the collection/extermination 
of toads during the colder months, as the vents act as 
a lure. Recent advances in cane toad control using 
fences (Florance et al., 2011) and pheromones (an 
alarm pheromone and an attractant pheromone can 
be used respectively to selectively kill or attract the 
tadpoles) (Crossland et al., 2012; Crossland & Shine, 
2012) may provide solutions. 


Further research and monitoring is required to better 
establish the threat the cane toad poses to GAB spring 
communities. This includes establishing an accurate 
distribution map of cane toads at GAB springs, as well 
as establishing if the cane toads are breeding within 
the springs. In addition to cane toads, the viability of 
endemic spring species is also at risk due to multiple 
additional threats, including mosquitofish (Gambusia 
holbrooki) and feral pigs (Fensham et al., 2010). 
Further research to inform the ongoing management 
of each individual threat is required to ensure 
conservation of these unique ecosystems. 


ACKNOWLEDGEMENTS 

This work was performed as part of the Queensland 
Government’s aquatic monitoring activities with 
permission from Bush Heritage Australia under the 
ethics permit DERM/201 1/05/02 and General Fisheries 
Permit 89474. Our thanks go to Adam Kerezsy from 
Bush Heritage Australia for his assistance with the 
collection of the cane toads and site access. 


LITERATURE CITED 

ALFORD, R.A., COHEN, M.P., CROSSLAND, 
M.R., HEARNDEN, M.N., JAMES, D. & 
SCHWARZKOPF, L. 1995. Population biology 
of Bufo marinus in northern Australia. Pp. 173- 
181. In Finlayson, M. (ed). Wetland Research in 
the Wet-Dry Tropics of Australia. (Office of the 
Supervising Scientist Report 101, Canberra). 

BLESSING, J.J., NEGUS, P.M. & CLIFFORD, S.E. 
2012. Spring Monitoring 2011: Great Artesian 
Basin Resource Operations Plan. (Queensland 
Department of Science, Information Technology, 
Innovation and the Arts, Brisbane). 

BUREAU OF METEOROLOGY 2012. Climate 
Statistics for Australian locations. (http://www. 
bom.gov.au/climate/averages/tables/cw_036007. 
shtml Accessed July 2012). 

CATLING, P.C., HERTOG, B.R.J., WOMBEY, J.C. & 
FORRESTER, R.I. 1999. The short-term effect of 


cane toads (Bufo marinus) on native fauna in the 
Gulf Country of the Northern Territory. Wildlife 
Research 26: 161-185. 

CROSSLAND, M.R. & ALFORD, R.A. 1998. 
Evaluation of the toxicity of eggs, hatchlings and 
tadpoles of the introduced toad Bufo marinus 
(Anura: Bufonidae) to native Australian aquatic 
predators. Australian Journal of Ecology 23: 
129-137. 

CROSSLAND, M.R., HARAMURA, T., SALIM, 
A.A., CAPON, RJ. & SHINE, R. 2012. 
Exploiting intraspecific competitive mechanisms 
to control invasive cane toads (Rhinella marina). 
Proceedings of the Royal Society B. 279(1742): 
3436-3442. 

CROSSLAND, M.R. & SHINE, R. 2012. Embryonic 
exposure to conspecific chemicals suppresses 
cane toad growth and survival. Biology Letters 8: 
226-229. 

DOODY, J.S., GREEN, B., RHIND, D., 
CASTELLANO, C.M., SIMS, R. & ROBINSON, 
T. 2009. Population-level declines in Australian 
predators caused by an invasive species. Animal 
Conservation 12: 46- 53. 

EHP. 2009. Climate change impacts in Queensland’s 
regions. The State of Queensland (Department 
of Environment and Heritage Protection). http:// 
www.ehp.qld.gov.au/climatechange/regional- 
summaries.html (accessed 10 April 2013). 

FAIRFAX, R, FENSHAM, R, WAGER, R, BROOKS, 
S, WEBB, A & UNMACK, P. 2007. Recovery of 
the red-finned blue-eye: an endangered fish from 
springs of the Great Artesian Basin. Wildlife 
Research 34: 156-166. 

FENSHAM, RJ. & FAIRFAX, RJ. 2009. 
Development and trial of a spring wetland 
monitoring methodology in the Great Artesian 
Basin. (Queensland Department of Environment 
and Resource Management, Brisbane). 

FENSHAM, R.J., PONDER, W.F. & FAIRFAX, 
R.J. 2010. Recovery plan for the community of 
native species dependent on natural discharge 
of groundwater from the Great Artesian Basin. 
Report to Department of the Environment, Water, 
Heritage and the Arts, Canberra. (Queensland 
Department of Environment and Resource 
Management, Brisbane). 

FENSHAM, R.J., SILCOCK, J.L., KEREZSY, A. & 
PONDER, W. 2011. Four desert waters: Setting 
arid zone wetland conservation priorities through 
understanding patterns of endemism. Biological 
Conservation 144 (10): 2459-2467. 


24 PROCEEDINGS OF THE ROYAL SOCIETY OF QUEENSLAND 


FLORANCE, D., WEBB, J.K., DEMPSTER, T., 
KEARNEY, M.R., WORTHING, A. & LETNIC, 
M. 2011. Excluding access to invasion hubs can 
contain the spread of an invasive vertebrate. 
Proceedings of the Royal Society B. 278(1720): 
2900-2908. 

FREELAND, W.J. 1984. Cane toads: a review of 
their biology and impact on Australia. (Parks and 
Wildlife Unit, Conservation Commission of the 
Northern Territory, Winnellie). 

FREELAND, W.J., DELVINQUEIR, B.L.J. & 
BONNIN, B. 1986. Food and parasitism of the 
Cane Toad, Bufo marinus, in relation to time since 
colonization. Australian Wildlife Research 13(3): 
489-499, 

GREENLEES, M.J., BROWN, G.P., WEBB, J.K., 
PHILLIPS, B.L. & SHINE, R. 2006. Effects of 
an invasive anuran (the cane toad, Bufo marinus) 
on the invertebrate fauna of a tropical Australian 
floodplain. Animal Conservation 9: 431-438. 

HEISE-PAVLOV, S.R. & LONGWAY, L.J. 2011. Diet 
and dietary selectivity of Cane Toads (Rhinella 
marina) in restoration sites: a case study in 
Far North Queensland, Australia. Ecological 
Management & Restoration 12: 230-233. 

HINKLEY, A.D. 1963. Diet of the Giant Toad, Bufo 
marinus (L.), in Fiji. Herpetologica 18: 253-259. 

HYSLOP, E.J. 1980. Stomach contents analysis - a 
review of methods and their application. Journal 
of Fish Biology 17: 411-429. 

KEARNEY, M., PHILLIPS, B.L., TRACY, C.R., 
CHRISTIAN, K.A., BETTS, G. & PORTER, 
W.P. 2008. Modelling species distributions 
without using species distributions: the cane toad 
in Australia under current and future climates. 
Ecography 31: 423-434. 

KEREZSY, A. 2011. Edgbaston Reserve: Recovery of 
red-finned blue-eye: Progress and recommendations. 
(Bush Heritage Australia Report). 

LEVER, C. 2001. The cane toad: the history and 
ecology of a successful colonist. (Westbury 
Academic and Scientific Publishing, West 
Yorkshire, United Kingdom). 

LINNAEUS, L. 1758. Systema naturae per regna 
tria naturae, secundum classes, ordines, genera, 
Species, cum _ characteribus,  differentiis, 
synonymis, locis. Tomus I. Editio decima, 
reformata. Holmiae. (Laurentii Salvi). 

PHILLIPS, B.L., BROWN, G.P. & SHINE, R. 2003. 
Assessing the potential impact of cane toads 
Bufo marinus on Australian snakes. Conservation 
Biology 17: 1738-1747. 


PONDER, W.F. & CLARKE, G.A. 1990. A radiation 
of hydrobiid snails in threatened artesian springs 
in Western Queensland. Records of the Australian 
Museum 42: 301-363. 

PONDER, WF, VIAL, M.,  JEFFERYS, 
E. & BEECHEY, D. 2010. The aquatic 
macroinvertebrates in the springs on Edgbaston 
Station, Queensland: A report for Bush Heritage 
Australia. (Australian Museum, Sydney). 

REED, R.N., BAKKEGARD, K.A., DESY, G.E. & 
PLENTOVICH, S.M. 2007. Diet composition of 
the invasive cane toad (Chaunus marinus) on Rota, 
Northern Mariana Islands. Pacific Conservation 
Biology 13(3): 219-222. 

SHINE, R. 2010. The ecological impact of invasive 
cane toads (Bufo marinus) in Australia. Quarterly 
Review of Biology 85: 253-291. 

SLUYS, R., GRANT, LJ. & BLAIR, D. 2007. 
Freshwater planarians from artesian springs 
in Queensland, Australia (Platyhelminthes, 
Tricladida, Paludicola). Contributions to Zoology 
76(1): 9-19. 

SOMAWEERA, R. & SHINE, R. 2012. The (non) 
impact of invasive cane toads on freshwater 
crocodiles at Lake Argyle in tropical Australia. 
Animal Conservation 15:152-163. 

SPEARE, R. 1990. A review of the diseases of the 
cane toad, Bufo marinus, with comments on 
biological control. Australian Wildlife Research 
17: 387-410. 

STRUSSMANN, C., RIBEIRO DO VALE, MB., 
MENEGHINI, M.H. & MAGNUSSON, 
W.E. 1984. Diet and Foraging Mode of Bufo 
marinus and Leptodactylus ocellatus. Journal of 
Herpetology, 18(2): 138-146. 

SUTHERST, R.W., FLOYD, R.B. & MAYWALD, 
GF. 1995. The Potential Geographical 
Distribution of the Cane Toad, Bufo marinus L. 
in Australia, Conservation Biology 10(1): 294- 
299. 

URBAN, M., PHILLIPS, B.L., SKELLY, D.K. & 
SHINE, R. 2007. The cane toad’s (Chaunus marinus) 
increasing ability to invade Australia is revealed by 
a dynamically updated range model. Proceedings of 
the Royal Society B 274: 1413-1419. 

WAGER, R. & UNMACK, PJ. 2000. Fishes of 
the Lake Eyre Catchment of central Australia. 
(Queensland Fisheries Service. Queensland 
Department of Primary Industries, Brisbane). 

WAGER, RN.E. & UNMACK, PJ. 2004. 
Threatened fishes of the world: Scaturiginichthys 
vermeilipinnis (Ivantsoff, Unmack, Saeed & 


DO CANE TOADS (RHINELLA MARINA) IMPACT DESERT SPRING ECOSYSTEMS? 25 


Crowley 1991) (Pseudomugilidae). Environmental Bufo marinus: a natural history resumé of 
Biology of Fishes 70: 330. native populations. Smithsonian Contributions 
ZUG, G.R. & ZUG, P.B. 1979. The marine toad, to Zoology 284: 1-54. 
AUTHOR PROFILES 


The authors are aquatic ecologists with an interest in assessing threats to freshwater ecosystems. Sara Clifford, 
Alisha Steward, Joanna Blessing and Peter Negus (project leader) undertake the Stream and Estuary Assessment 
Program (SEAP) within the Water Planning Ecology Group (WPE) of the Queensland Department of Science, 
Information Technology, Innovation and the Arts. SEAP is a state-wide program which monitors the condition 
and trend of riverine ecosystems, specifically relating to anthropogenic risks to condition whilst confirming 
and improving on current understanding of processes. As Principal Scientist of WPE, Jon Marshall is actively 
involved in all of the group’s projects, including SEAP. 


STYGOFAUNA PRESENCE WITHIN FRESH AND HIGHLY SALINE AQUIFERS 


OF THE BORDER RIVERS REGION IN SOUTHERN QUEENSLAND 
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Stygofauna are minute, inconspicuous biota that inhabit aquifers and have remained largely unsampled 
in Queensland until recently. The presence of stygofauna in some aquifers has been considered unlikely 
due to extreme physico-chemical characteristics such as high electrical conductivity. However, until 
explored, the presence or absence of stygofauna in such aquifers cannot be confirmed. 


Stygofauna were sampled from fifteen monitoring bores using phreatobiological nets in six sub-catch- 
ments of the Border Rivers region in southern Queensland. Stygofauna from three taxa (Bathynellacea, 
Cyclopoida, Harpacticoida) were found at most sites, from a variety of groundwater systems. Ad- 
ditionally, stygofauna were collected from aquifers where they were not expected to occur, including 
three bores located in clay-dominated groundwater systems with very high electrical conductivities, 
approaching that of seawater, between 36,300 and 54,800 S/cm. 


The findings from this investigation show that the presence of stygofauna is not always limited by 
the physico-chemical parameters of aquifers, particularly electrical conductivity. Queensland contains 
extensive groundwater resources that are largely unexplored in terms of stygofauna, and it is likely that 
future investigations into these aquifers will reveal rich and diverse groundwater ecosystems. 
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INTRODUCTION 

Ninety-seven per cent of the world’s unfrozen 
freshwater is found within subterranean aquifers 
(Gibert et al., 1994). Traditional thinking considered 
aquifers as vast, lifeless, storage areas (Hancock & 
Boulton 2008), but this view is changing as ecologists 
and microbiologists intensify their exploration of these 
environmental systems (Gibert & Deharveng, 2002; 
Boulton et al., 2008). The interstitial spaces within 
aquifers can be habitats for bacteria, as well as obligate 
groundwater invertebrates known as “stygofauna” 
(Tomlinson et al., 2007). Tomlinson et al. (2007) 
proposed a number of values associated with stygofauna 
communities: as a reserve of biodiversity; as potential 
providers of ecosystem goods and services; and as 
biological indicators of the “health” of groundwater 
ecosystems. Bacteria and stygofauna may also play an 
important role in groundwater purification and carbon 
and energy cycles (Serov, 2002). Bacteria are used in 
the removal of toxins from contaminated groundwater 
(‘bioremediation’, Chapelle, 1999), while stygofauna 
can alter the interstitial pore size through their 
movement, maintaining hydraulic flow pathways and 
transporting material through groundwater ecosystems 
(Hancock et al., 2005; Creuze et al., 2009). 


The extreme conditions imposed on _ stygofauna 
through an entirely subterranean existence are 
believed to place strict boundaries on_ their 
distribution in groundwater systems (Hancock et 
al., 2005). Dole-Olivier et al. (2009) indicated that 
stygofauna biodiversity hotspots mainly occurred 
in highly-permeable geological formations such as 
karst and coarse alluvial environments. The highly 
porous substrate of these environments allows for 
the movement of individuals within the aquifer and 
facilitates the movement of high quality recharge 
water delivering oxygen, carbon and nutrients to the 
aquifer and the ecosystem (Galassi et al., 2009). 


There has been significant literature written on 
general landscape drivers and the attributes of 
habitats supporting subterranean ecosystems 
(Hancock et al., 2005, Humphreys, 2006; Boulton 
et al., 2008). Despite broad consistencies observed 
across stygofauna communities, numerous studies 
have failed to identify specific physical or chemical 
groundwater quality variables that demonstrate clear 
relationships with stygofauna diversity or abundance 
(Boulton et al., 2008). Hancock & Boulton (2008) 
found greater taxa richness in east Australian sand 
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and gravel alluvial aquifers associated with shallow 
water tables and electrical conductivities less than 
1500 uS/cm. Proximity to deep-rooted phreatophytic 
vegetation has also been suggested as a valuable 
source of organic carbon to subterranean ecosystems 
(Hancock & Boulton, 2008; Dole-Olivier et al., 2009). 
These habitat attributes are believed to influence the 
diversity and abundance of stygofauna populations. 


Australian groundwater systems were originally 
thought to be deficient in subterranean fauna compared 
to other regions of the world due to a relative paucity 
of karstic habitats and widespread aridity (Guzik et 
al., 2010). However, as karst and cave systems were 
explored, increasingly rich and diverse subterranean 
ecosystems were uncovered (Humphreys, 1999; 
Humphreys, 2006; Watts & Humphreys, 2006). As the 
higher than expected subterranean biodiversity became 
evident, Australia has more recently been predicted to 
contain a wealth of stygofauna populations due to the 
exceptional diversity of groundwater systems present 
(Guzik et al., 2010). 


The majority of stygofauna sampling in Australia has 
been conducted in the arid western half of the continent 
(Guzik et al., 2010). The concentration of research 
effort originated due to the presence of karst systems 
in Western Australia. Once significant subterranean 
diversity was established, survey effort was reinforced 
through environmental regulations _—srequiring 
subterranean biological surveys associated with the 
environmental review process for the mining industry 
(Guzik et al., 2010, Hancock & Boulton, 2008). 


Stygofauna populations have since been found in 
most Australian states and in a diverse array of 
groundwater systems including alluvium, anchialine 
pools, calcretes, fractured rock, karst, pseudokarst, 
and pisolites (Guzik et al., 2010). The Western 
and South Australian governments have both 
conducted extensive biological surveys to document 
stygofauna diversity. Stygofauna research in eastern 
Australia is developing through some PhD projects 
(e.g. Hancock, 2004a; Tomlinson, 2008) and the 
recent requirement for stygofauna investigations 
during the Environmental Impact Statement (EIS) 
process. There is growing acknowledgment of the 
importance of Groundwater Dependant Ecosystems 
(GDEs) including stygofauna across Queensland, 
and the need to improve the collation of knowledge 
and understanding (Tomlinson & Boulton 2008; 
Tomlinson, 2009; Perna, et al., 2012). 


In Queensland, stygofauna distributions are known 
from a limited number of surveys, mainly in coastal 
catchments (Hancock & Steward, 2004; Hancock, 
2004a; Hancock, 2004b; Tomlinson, 2008), and some 
opportunistic sampling within irrigated agricultural 
land west of the Great Dividing Range (Tomlinson, 
2009). Although investigations to date remain 
relatively limited in the spatial coverage across 
Queensland groundwater systems, a rich diversity of 
stygofauna is starting to be uncovered. As the mining 
and gas sectors develop in Queensland, particularly 
in the semi-arid central western catchments, there is 
increasing need for deficiencies in knowledge and 
understanding of stygofauna to be addressed. 


As our awareness of Australian subterranean 
ecosystems develops, there is a growing appreciation 
of the biological assets still to be documented (Hancock 
& Boulton 2008). The high number of groundwater 
systems in Queensland yet to be investigated suggests 
further discoveries are likely as exploration continues. 
General patterns in research to date can guide the 
exploration of new areas; however, the prevalence 
of exceptions within this developing field advocates 
regular testing of assumptions concerning patterns 
in habitat attributes. The current study aims to 
confirm the presence of stygofauna populations in 
an unstudied region, and collect baseline information 
on population distributions to guide future targeted 
biological surveys. 


MATERIALS AND METHODS 

STUDY AREA AND CATCHMENT CONDITIONS 

The Border Rivers catchment in the upper Murray 
Darling Basin, Australia, is located on either side of 
the State border between Queensland and New South 
Wales, immediately west of the Great Dividing Range. 
The catchment covers an area of 63,550 km/’, and 
lies in approximately equal portions in northern New 
South Wales and southern Queensland. The border 
passes along the parts of the Barwon, Macintyre and 
Dumaresq Rivers, with the Queensland Border Rivers 
catchments including the Severn River, Pike Creek, 
Macintyre River, Macintyre Brook and the Weir River. 
The groundwater systems within the study area range 
from highly productive alluvial sands and gravels, 
fractured granite, traprock alluvium and highly saline 
clay-dominated systems (Figure 1). 


The region often experiences periods of prolonged 
drought followed by extreme flooding. In the six months 
preceding the current study, the Border Rivers catchments 


STYGOFAUNA PRESENCE WITHIN FRESH AND HIGHLY SALINE AQUIFERS 29 
OF THE BORDER RIVERS REGION IN SOUTHERN QUEENSLAND 


416200264 4162700154 


~~ ne, 


O Rivet 416200044 
Pe i Cer 


@ Border River Fractured Granite 
@ Sracker Creek Alluvium 

@ Commoron Creek Alluvium 

O Dumaresg Alluvium 

@ Macintyre Brook Alluvium 

@ Macintyre River Alluvium 

© Traprock Alluvium 

© weir Atuvium 

o 0 _ 





FIGURE 1. The distribution of sampling sites in the Border Rivers catchment, defined by groundwater system. 
Note: sites in close proximity may overlap on map. Inset — the location of the Border Rivers catchment in Queensland. 


experienced some of the highest rainfalls and flooding 
on record (Figure 2). The influence on stygofaunal 
populations of such excessive and prolonged flooding is 
difficult to determine. In addition to significant surface 
water flows, it provides a considerable groundwater 
recharge event across the region. Furthermore, there is 
limited understanding of the hydraulic and ecological 
time lags associated with aquifer ecosystem responses to 
events of this magnitude. 


SITE SELECTION 

Hancock & Boulton (2008) investigated stygofauna 
distribution in Eastern Australian alluvial aquifers, 
suggesting abundance and diversity were related to 
the efficiency of recharge water to deliver oxygen and 
nutrients to the ecosystem, the size of void spacing 
and less extreme water quality conditions. 


Aquifers were selected for sampling along a gradient of 
preferred habitat characteristics based on the currently 
available literature, and included the full diversity of 


groundwater systems present, even those considered 
as less favourable for subterranean fauna (Hancock 
& Boulton 2008). This approach to site selection 
was intended to collect information on patterns of 
stygofauna population distribution across the region. 


Significant consultation with local groundwater experts 
targeted the selection of monitoring bores across 
all groundwater systems and the diversity of bore 
characteristics. Sites were selected along a gradient 
of habitat characteristics based on depth to water 
table, hydraulic flow rates, void spacing (substrate 
types) and water chemistry. The resulting sampling 
network contained fifteen Queensland Government 
groundwater monitoring bores distributed to represent 
the diversity of groundwater systems found within the 
Border Rivers study area (Figure 1). 


SAMPLING TECHNIQUE 
Stygofauna sampling was conducted over two 
weeks in May 2011, the beginning of the dry season. 
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FIGURE 2. River flow plot from a representative gauging 
station within the study area (GS 416201A, Macintyre River 
at Goondiwindi). Stygofauna sampling occurred during May 
2011. This flow plot is included to demonstrate the extreme 
rainfall experienced across the region in the six months 
preceding the stygofauna investigation. 


Stygofauna samples were collected using a weighted 
phreatobiological net with 50 um mesh and 50 mm 
diameter. Nets were lowered to the bottom of a bore 
and then drawn 50 cm through the water column four 
times to agitate resting fauna before a slow retrieval. 
Net contents were transferred to a 50 um sieve. 
Sampling was repeated six times as per Hancock & 
Boulton (2008) at each bore, with collected material 
transferred into vials, stored in 100 % ethanol and 
stained with Rose Bengal dye (as demonstrated 
in Figure 3) to aid in the detection of fauna during 
laboratory processing. 


Groundwater quality parameters were sampled at 12 of 
the 15 bores where pumping was logistically possible. 
Bores were pumped using a “Waterra power-pump 2’, 
removing up to 300 L of groundwater into 10 L containers. 
Groundwater quality was monitored from the final 10 L 
collected, with hand held water quality probes for electrical 
conductivity, dissolved oxygen and pH, and parameter 
values recorded after having stabilised. Care was taken to 
ensure water exited the pump below the water level of the 
bucket so that aeration was minimised. Three bores were 
unable to be pumped due to time restrictions and poor 
pumping rates— from these bores water quality information 
was sourced from the most recent routine monitoring 
data stored in the Queensland Government groundwater 
database. Groundwater system characteristics for each 
bore were sourced from drillers’ records and bore cards 
stored in the Departmental database. 





FIGURE 3. Photograph of Bathynellacea specimens. The 
specimens are stained pink using Rose Bengal dye to assist 
in sample processing, as stygofauna typically lack pigment. 


The monitoring effectiveness of stygofauna net 
sampling compared to pump sampling of bores has 
been discussed in the literature (Hancock & Steward, 
2004; Hancock & Boulton, 2009). Hancock & 
Boulton (2009) indicated pump sampling collected 
more taxa and higher total abundances than net 
sampling. However, net sampling is advantageous 
under specific circumstances due to lower costs, 
faster sampling and applicability to a wider range of 
bores, such as low yield groundwater systems. The 
preliminary nature of the current study supported the 
use of net sampling equipment, due to the logistic 
advantages of the greater number of sites, range of 
bores and groundwater systems able to be sampled. 


STYGOFAUNA 
ANALYSIS 

Stygofauna sample contents were rinsed to remove 
excess Rose Bengal stain, and then sorted under a 
stereo-microscope (126 x magnification). Specimens 
were identified to the lowest taxonomic level possible, 
and enumerated. Often, order or family level is the 
only taxonomic resolution achievable for stygofauna, 
as many Australian taxa are new to science and remain 
undescribed. For the purposes of this paper, only 
records of Copepoda (Harpacticoida and Cyclopoida) 
and Syncarida (Bathynellacea) are reported. We 
could not confidently identify other groups, such as 
Nematoda and Acarina, as being ‘true’ stygofauna 
taxa, as they can also occur in surface waters and 
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FIGURE 4. The proportion of each stygofauna taxon, based on abundance data, found in each sample in the Border Rivers catchment. 


may have dispersed to the sampled groundwater 
systems during the extreme flooding event preceding 
sampling. The proportion of each stygofauna taxon in 
each sample was calculated based on abundance data, 
and plotted as a pie graph. 


RESULTS 

Stygofauna were present at multiple locations across 
multiple groundwater systems in the Border Rivers 
catchment. In total, 694 stygofaunal invertebrates 
across 3 taxonomic orders (Bathynellacea, Cyclopoida, 
and Harpacticoida) were collected from bores during 
this investigation (Table 1). The spatial distribution of 
stygofauna results is displayed in Figure 4. The number 
and diversity of taxa found in any one bore varied 
greatly, with some bores limited to a single individual 
and one bore in the Dumaresq River alluvium containing 
579 individuals across 2 orders. One bore contained 
individuals from all three taxonomic orders. Six bores, 
including sites in the Macintyre Brook, Macintyre 
River, Border River, Weir River, and Dumaresq River, 
did not contain stygofauna. 


Bore and groundwater system attributes varied across 
the sampled bores (Table 2). The depth to water table 
ranged between -0.75 m and -21.30 m, with stygofauna 
present at depths between -0.85 m and -21.30 m. 
The electrical conductivity of sampled groundwater 
ranged between 150 uS/cm and 54,800 uS/cm, with 
stygofauna present at conductivities between 309 uS/ 
cm and 54,800 uS/cm. pH ranged between 3.5 and 8.1, 
with stygofauna present in pHs between 3.5 and 7. The 
sampled groundwater system substrates varied greatly, 
ranging from highly productive alluvial sands and 
gravels, fractured granite, traprock alluvium and highly 
saline clay-dominated systems. Stygofauna were found 
to be present across the range of substrate types, except 
for the Macintyre Brook and Macintyre River alluvia. 


DISCUSSION 
This investigation into subterranean fauna in the 
Border Rivers confirmed the presence of stygofauna 
populations in multiple groundwater systems across 
the Border Rivers catchments. There was no obvious 
Spatial pattern in the presence or absence of stygofauna, 
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TABLE 1. Abundance of stygofauna collected from bores in the Border Rivers catchment. 
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EC 
ESC 
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ramaren [Masi ookaTinwm YO 





with Bathynellacea, Cyclopoida and Harpacticoida 
found throughout the catchment. An absence of 
stygofauna from aquifer samples is not uncommon 
in the literature (e.g. Deharveng et al., 2009; Hahn 
& Fuchs, 2009; Martin et al., 2009). All three orders 
were present in the upper and lower parts of the 
catchment, and Bathynellacea and Cyclopoida were 
found in both fractured rock and alluvium aquifers. 
Such heterogeneous distributions have been seen in 
European groundwater systems, with ‘hotspots’ of 
taxa found in some areas, and no stygofauna found 
in others (Deharveng et al., 2009; Martin et al., 
2009). The diversity of groundwater systems found 
to support stygofauna populations highlights the need 
for future investigations to consider all aquifer types 
in any quantitative survey of abundance or diversity. 


The presence of stygofauna in the highly-saline, clay- 
dominated Weir River alluvium, with conductivities 
between 39,000 pS/cm and 54,800 uS/cm, was an 
unexpected result. It has been suggested that inland 
stygofauna have a preference for lower salinities, 
with higher abundances and diversity found at 


conductivities less than 1500 uS/cm (Hancock & 
Boulton, 2008). However, stygofauna from all 
three taxonomic orders were found in systems 
approaching or exceeding marine conductivities, 
suggesting wider tolerances for physical conditions 
than previously expected. 


Additionally, the high clay content of the Weir 
River alluvium and the associated low hydraulic 
conductivities and lack of interstitial space was not 
expected to provide habitat to support stygofaunal 
communities. It has been suggested in the literature 
that phreatophytic tree roots may constitute a unique 
habitat, creating larger interstitial spaces and providing 
a source of organic material for subterranean fauna 
(Hancock & Boulton, 2008). The Weir River alluvium 
sites were located alongside areas of relatively 
natural land use such as grazing or re-vegetated road 
reserves, with close proximity to relatively dense tree 
populations. The influence of phreatophytic tree roots 
was not specifically investigated in the current study, 
but is suggested as a likely contributing factor for 
stygofauna presence in such extreme habitats. 
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The use of net sampling in the Border Rivers 
has provided a valuable qualitative indication of 
stygofauna presence and distribution. Net samples 
alone have been found to under-sample groundwater 
faunal communities, with comprehensive assessments 
of biodiversity requiring a combination of net and 
pump samples (Hancock & Boulton, 2009). Repeated 
sampling in the Border Rivers with the addition of 
pump sampling techniques is expected to uncover 
additional taxa and improved abundance estimates. 
Future quantitative biological surveys will describe 
the richness of Border Rivers stygofauna communities 
compared to other groundwater systems across 
Australia. 


During the six months preceding sampling, the 
Border Rivers experienced an extreme flooding 
event. Anticipating the influence of such intense 
hydraulic conditions is difficult, particularly without 
comparable pre-flood population data. High flow 
events are known to recharge groundwater systems, 
and often provide delivery of oxygen, nutrients and 
carbon essential to their ecosystems. Time lags 
associated with recharge delivery is influenced by 
substrate characteristics and hydraulic conductivity 
traits. Additionally, the population dynamics of low 
metabolic stygofaunal ecosystems and the time taken 
for energy and nutrient inflows to be converted into 
population abundance is unknown. Such an extreme 
flood event is also likely to have over-topped some 
bore casings. The impact of such direct inflows of 
surface waters to aquifer systems, and how this may be 
interpreted through measured parameters is unknown 
(Hancock et al., 2005). The current study does not 
attempt to predict or interpret such influences, but 
notes their occurrence for comparison with results 
from future investigations. Repeated sampling to 
address temporary variability is recommended (Dole- 
Olivier et al., 2009), and has yielded additional taxa in 
other studies (Hancock & Boulton, 2009). 


CONCLUSIONS 

This investigation highlights stygofauna as a diverse and 
adaptable taxonomic group that is able to occupy a broad 
range of environmental conditions including aquifers 
with very high salinity. Future monitoring of stygofauna 
should not be limited by untested, pre-existing ideas of 
distribution or habitat preferences. Queensland contains 
extensive groundwater resources that are largely 
unexplored in terms of stygofauna. It is likely that future 
investigations into these aquifers will reveal them to be 
rich and diverse groundwater ecosystems. 
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BASIN, USING STATISTICAL AND CONCEPTUAL METHODS 


McNEIL, V.H. & RAYMOND, M.A.A. 


In Australia the precautionary principle is one of the guiding tenets of ecologically sustainable develop- 
ment, and implies the need to define acceptable limits for water chemistry, if necessary on the basis of 
current scientific information. Groundwater chemistry is difficult to characterise with certainty because 
of the many interacting controls, which may contribute to the typically high variability. A desktop 
methodology is presented, which partitions a large, diverse catchment, the Fitzroy River Basin, into 
reasonably homogeneous groundwater chemistry zones on the basis of an extensive historical dataset, 
collected opportunistically at irregular intervals. A dual approach was adopted for assessment, with 
multivariate analysis supported conceptually by a range of environmental factors and presented as GIS 
background layers. Surface water samples were included in the multivariate analysis to define areas 
of interaction. Assessment was based on salinity and the major ions, which are reasonably insensitive 
to sampling and storage protocols, and the data were checked for reliability through chemical balance 
and independent salinity measures. The data were sorted into water types using cluster and principal 
component analysis, and then bores were classified with their predominant water type and plotted as 
a map layer for visual comparison with conceptual criteria of geology, climate and land use. Spatial 
bias in the sampling necessitated subjective judgement in delineation of zonal boundaries, decisions 
on degree of subdivision, and explanation of uncertainty across the basin. This enabled the definition 
of 44 discrete groundwater chemistry zones in the Fitzroy Basin. Although the study is too broad scale 
to evaluate processes, it indicated that the groundwater contains two major chemical sequences, each 
evolving through a wide range of salinity. One is consistent with chemical weathering of Palaeozoic or 
basaltic rocks or derived alluvium in humid areas, and the other is a more sodic sequence associated 
with sub-coastal, quartz rich terrains such as sandstones and granites. 


Vivienne H. McNeil (Vivienne. McNeil@qld.gov.au), and Myriam A.A. Raymond (Myriam.Raymond@ 
qld.gov.au) Water Planning Ecology, Science Delivery, Queensland Department of Science, Informa- 
tion Technology, Innovation and the Arts, GPO Box 2454, Brisbane OLD 4011, Australia. 
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MAPPING REGIONAL GROUNDWATER CHEMISTRY ZONES IN THE FITZROY 


INTRODUCTION 

Salinity related water quality problems, including the 
composition of individual salts, can seriously affect 
agricultural production, the viability of infrastructure, 
and the welfare of regional communities (NAP 2002). 
They are also a priority issue for maintaining ecological 
health and biodiversity (Dunlop & McGregor 2007). 
Groundwater is a major and increasingly significant 
resource in Queensland, particularly in rural areas, 
and supports a range of groundwater dependent 
ecosystems. It also interacts with soils and surface 
waters, either naturally or as a result of human 
activities. Despite its importance, there is limited 
knowledge and understanding of the resource over 
wide areas of Queensland because of the size of the 
state, and its low but increasing population density. 


The precautionary principle is one of the guiding tenets 
of ecologically sustainable development in Australia, 
and is included in a variety of legislation and policies 
(Peterson (2006). In order to ensure groundwater 


protection under this principle, acceptable or ‘baseline’ 
limits for groundwater chemistry must be established 
on the basis of current scientific information to 
support adaptive management pending more intensive 
investigations where necessary. Baseline limits have 
been defined by Shand et al. (2007) and Edmunds 
et al. (2003) as the range in concentration which 
occurs through natural processes within a particular 
groundwater system. Much of the current scientific 
information on Queensland groundwater chemistry 
consists of large, but unevenly collected and poorly 
documented analysed samples. 


Salinity and water chemistry zones with baseline 
ranges have previously been established for 
Queensland’s _surface-waters, using _ statistical 
techniques designed for bulk processing of historical 
data (DERM 2009), but defining — sufficiently 
homogeneous chemical zones in groundwater is a 
more complex task because water sources are less 
clear, spatial variation is three dimensional, and many 
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factors influence the chemistry. Such factors include 
recharge composition, soils, geology and rainfall, as 
well as the localised influence of human activities and 
the interactions between water bodies. The resulting 
high natural variability of groundwater chemistry 
may breach limits for environmental values, even in 
the absence of human impacts (Shand et al. 2007; 
McNeil, Cox & Preda 2005; Hancock 2002; Searl, 
Armitage & Dawson 1999). 


The task of spatially partitioning groundwater has 
been frequently addressed in the literature. ANRA 
(2000), for instance, defined an Australia-wide 
system of Groundwater Management Units (GMUs) 
which refer to hydraulically connected regional 
groundwater systems. Victoria has further defined 
groundwater provinces, and produces a detailed guide 
of condition and vulnerability (ANRA 2002). Some 
large scale projects have applied more intensive 
methods to determine baseline water quality ranges, 
for instance Shand et al. (2007) and Eberts & George 
(2000) used datasets which included trace elements 
and stable isotopes, and compiled detailed reviews 
of local catchment information, but this is beyond 
the resources available for many studies of larger 
catchments. Most of these authors acknowledge that 
multivariate analysis of water quality data should 
be supplemented by conceptual information when 
partitioning groundwater bodies. To this end, a number 
of computational studies in self-organised mapping 
have been carried out in a variety of research areas 
where more rapid spatial mapping is required, for 
example Corchado & Perez (2011), Csaky & Please 
(2003), Di Prinzio, Castellarin & Toth (2011), Whelan 
et al. (2010) and Lamble & Fraser (2005). 


Changes in water chemistry baselines with depth, 
as discussed by Acworth & Janowski (1993), is 
well recognised in regional groundwater studies, 
for instance Eberts & George (2000); Shand et al. 
(2007). Salinity in undisturbed aquifers is usually 
reported to increase with depth below the water table, 
as water residence time increases and rock/water 
reactions progress (McMahon, Bohlke & Lehman 
2004, Kumar & Sinha 2010), although this is not 
always the case. For instance, McMahon, Bohlke & 
Lehman (2004) found a band of higher salinity close 
to the surface of the water table of the Texas (USA) 
High Plains aquifer, a condition that they attributed 
to overlying agriculture. Similar instances have been 
reported by other authors such as Lambrakis (2006) 
and Rosenthal et al. (1992). Alternatively, in places 


where there are multiple discrete aquifers as in 
much of the Great Artesian Basin, individual aquifer 
chemistries may be unrelated to each other. 


Another concept relevant to the maintenance 
of chemical baselines is ground/surface water 
interaction, either through recharge or baseflow. 
Areas of interaction may be vulnerable to human 
activity through land clearing, transfer of pollutants, 
contamination by effluents, or changes to natural 
hydrological regimes (for example: McMahon, 
Bohlke & Lehman 2004; Karro, Marandi & Vaikm 
2004; and Rosenthal et al. 1992). 


This current study has combined and adapted 
approaches from the foregoing literature to group the 
groundwater data from a large and complex catchment 
into chemical water types, to enable zones of similar 
groundwater chemistry to be rapidly defined, even 
if the data is unevenly distributed in space and time. 
Zones were defined on the basis of predominant 
water types, in combination with environmental 
information provided through multiple GIS layers. 
Current baseline water chemistry was then estimated 
for each zone, based on the percentiles of its internal 
ranges. The addition of surface water samples to the 
groundwater data was introduced in order to identify 
areas of similarity and therefore possible interaction. 
Since no data is available pre-development, it is 
acknowledged that the zones and their baseline 
ranges only represent current mid-range levels, but 
these interim values are in line with the precautionary 
principle in providing a filter to identify outlying sites 
and sudden or rapid change. Groundwater models and 
more intensive assessment methods, such as those 
presented in Raiber, et. al. (in press) and Daughney 
et. al (2012), can be applied at a later date to identify 
contributory processes and refine ranges in selected 
localities of high priority. 


MATERIALS AND METHODS 

STUDY AREA 

The Fitzroy River Basin is a major inland river system 
in Central Queensland (Figure 1) with an area of 
approximately 143,000 km”. It is the largest catchment 
draining to the east coast of Australia, discharging into 
the southern end of the Great Barrier Reef Lagoon. 
There are five major sub-catchments, the Nogoa, 
Comet, Dawson, Isaac and Mackenzie, divided by a 
series of low mountain ranges and drained by streams 
which are deeply incised into their floodplains (Douglas 
et al. 2006). The Nogoa River rises in the far west, and 
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FIGURE 1: The Fitzroy River Basin with tributary sub-basins. 
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flows through Emerald to meet Theresa Creek from 
the north-west near Clermont. The Comet River drains 
the south-central catchment, and forms the Mackenzie 
River at its junction with the Nogoa. The Connors/Isaac 
system, which drains the northern part of the catchment, 
joins the Mackenzie just upstream of Tartrus Weir. This 
is then joined further downstream by the largest sub- 
catchment, the Dawson, which drains the southern part 
of the basin with the Don and Dee Rivers as its major 
tributaries (Douglas et al. 2006, Dutvenvoorden et al. 
2001). The coastal alluvium deposited by the Fitzroy 
River, 1s constrained to the north by steep hills of 
rhyolite and granite, including Mt Archer (608m ASL) 
Just outside Rockhampton, and the Berserker and Flat 
Top Ranges running south-east towards the estuary of 
the Fitzroy River, which has more extensive plains and 
wetlands on its southern side. The population of the 
Basin is about 155,000 (OESR 2000). Rockhampton 
is the largest city, other major towns being Emerald, 
Blackwater, Taroom, Theodore, Clermont, Wandoan, 
Injune and Biloela. Principal land uses include grazing, 
irrigation, mining, and power generation, and there are 
also significant areas of state forest and national park. 


The basin is characterised by a sub-tropical, semi- 
arid climate with high rainfall variability (Calvert et 
al. 2000), and a variety of groundwater systems are 
present from both alluvial and non-alluvial aquifers. 
There are also broad regions where groundwater is only 
found sporadically. Groundwater is regularly obtained 
from the alluvial flood plains of the Fitzroy and its 
tributaries, particularly the Dawson and Isaac sub- 
catchments, and from scattered bores in alluvium of the 
Comet River, the lower reaches of the Nogoa (mainly 
irrigated by surface water from Fairbairn Dam), and 
around Sandy Creek near Clermont (ANRA 2000). 
Most alluvial aquifers are shallow and unconfined, with 
average bore depths of 20 m or less. The waters tend 
to be hard and relatively saline, but usually suitable for 
stock and domestic use (McNeil & Raymond 2011). 
The most economically important alluvial aquifers are 
found in the Dawson sub-catchment, both upstream 
on the Dawson around Theodore, and on the Don and 
Dee tributary system, particularly in the Callide Creek 
area that is situated around Biloela behind the Coastal 
Range. These aquifers support the irrigation of cotton 
and lucerne (ANRA 2000). 


Other types of aquifers, shown in Figure 1, are also 
widespread in the Fitzroy, and have a role in supporting 
the pastoral industry, dairying, small crops and town 
supplies (Pearce 1982). The most important non- 


alluvial aquifers are two basalt flow systems, one 
extending from Emerald to Clermont in the north, and 
the other around Springsure in the south. The Springsure 
basalts are vesicular, and higher yielding than the 
northern flows where supplies are mostly obtained 
from fracture zones. In addition to the basalts, there are 
karstic aquifers in Devonian limestone lenses around 
Mt Larcom, and water yielding vesicular Permian 
andesites near Mt Morgan (Pearce 1982). Porous 
sandstones cover large areas of the basin, sometimes 
associated with coal measures (Ishaq 1985; Cadman et 
al. 1998). Those in the central area are related to the 
underlying Bowen Basin, while those in the south, 
around the headwaters of the Dawson and Nogoa 
Rivers, are part of the Surat Basin and form part of the 
Great Artesian Basin Eastern Recharge area (Fensham 
et al. 2007; Cadman et al. 1998; Baxter 1992; Loch & 
Rolfe 2000). Limited supplies of groundwater can also 
be obtained from unconsolidated Cainozoc sediments 
which overlies much of the landscape in the northern 
and central part of the basin. 


DATA 

More than 13,000 subartesian water quality samples 
have been collected by the Queensland Department 
of Natural Resource Management (DNRM) from 
4,780 bores since the 1950s. Of these, 9,659 are 
from shallow bores, which are regarded for reasons 
discussed later to be those with a depth of less than 
30m. Shallow groundwaters are considered more 
likely to be unconfined and to interact with the 
surface. In addition, 5,037 surface water samples 
from the Queensland Government surface water 
database were used to compare local ground and 
surface water chemistry. These were accumulated 
through ambient monitoring of 110 gauging stations 
in the basin, as well as for specific projects by various 
organisations such as the Australian Department of 
Natural Resources, Mines and Energy (NRM&E), the 
Australian Centre for Tropical Freshwater Research 
(Congdon 1991; McNeil & Cox 2000) and the Fitzroy 
National Landcare Program (Noble et al. 1997). 


The available water quality parameters almost always 
included electrical conductivity (EC), usually also 
accompanied by pH, colour, turbidity, Si0,, hardness, 
alkalinity, Sodium Adsorption Ratio (SAR), Residual 
Alkali Hazard (RAH), Total Dissolved Ions (TDI), 
Total Dissolved Solids (TDS), Na, K, Ca, Mg, Fe, 
Mn, HCO,, CO,, Cl, F, NO,, and SO,, accompanied 
by groundwater Level. All of the data were analysed 
at the NATA registered Queensland Health Scientific 
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Services Laboratories, with high levels of internal 
quality control and a documented history of changes 
in analytical procedures. However, because the 
sampling was unplanned and opportunistic, its spatial 
and temporal distributions are uneven. This results 
in differing levels of confidence in defining the 
boundaries of zones across the basin, and therefore in 
the accuracy of the baseline ranges. 


It is acknowledged that standardised sampling and 
preservation protocols were not introduced until 
the early 1990s, and that long storage times may 
have reduced the representativeness of older data in 
terms of reflecting aquifer conditions. However, it is 
known that the bores were sampled after purging by 
accredited hydrologists and also that the EC and major 
ions, on which this study depends, are sufficiently 
resilient to lack of sample preservation and extended 
storage to statistically define broadscale ranges in a 
large volume of data (McNeil 2002). 


APPROACH TO ASSESSMENT 

The purpose of this assessment is to sort the 
groundwaters of the entire catchment into spatial 
zones, large enough to be clearly defined and 
managed, but within which the water quality baseline 
ranges are not too broad to be relevant. Because of 
the uneven data coverage as described above, a dual 
statistical and conceptual approach has been adopted 
for zone definition. First, the water samples were 
sorted into chemical types using robust multivariate 
statistical techniques. Then spatial zones were 
subjectively defined on the basis of regional water 
chemistry, supported by a range of environmental 
factors superimposed as GIS layers. 


Before commencing the multivariate analysis, a series 
of tests were applied to the data to exclude incomplete 
or chemically inconsistent analyses. These tests were 
based on the chemical balance of the major cations 
(Na, K, Ca, and Mg) and anions (HCO,, CO,, Cl, 
NO,, and SO,), as well as by comparing independent 
salinity measures. The balance test used was the 
formula for the percent charge balance error (YoCBE) 
recommended by Freeze & Cherry (1979) and 
Eaton et al. (1995). The equation for the %CBE 1s: 


, 100*(5\z*m, —}.z*m,) 


I) 


where z = the charge on an individual anion or cation 
(absolute value), mc = molality of a cation, and ma 
= molality of an anion. Three salinity parameters 
were included in most of the data. TDI, and TDS 
are calculated measures, and were recomputed using 
other constituents. The TDI and the independently 
measured EC were then correlated to locate samples 
where the ratio was outside normal ranges as presented 
in McNeil & Cox (2000). 


EXPLORATORY MULTIVARIATE ANALYSIS 
After incomplete or unreliable analyses had been 
removed, the surface and groundwater chemistry 
data were collated and subjected to an exploratory 
multivariate analysis. The aim of this type of analysis 
is not to test a hypothesis or establish statistical 
significance, but to reveal patterns in large or complex 
datasets. Accordingly, the results cannot independently 
justify why such patterns exist (Ramette 2007). The 
suite of methods chosen for an exploratory analysis 
must be capable of relating the data to factors of 
interest, and this requires subjective judgement as 
alternative methods tend to produce different results. 
The multivariate methodology used here is similar to 
those that were used by many other hydrogeological 
studies as cited below, and past experience, for 
instance McNeil, Cox & Preda (2005) and McNeil 
(2002), has shown it to be a sound basis for classifying 
the bulk chemistry of Queensland natural waters in 
terms underlying geology and climatic patterns. 


As described in McNeil, Cox & Preda (2005), cluster 
analysis was selected as the primary assessment 
method as it is nonparametric and is suitable for 
the known distribution types within the data. The 
parameter ranges were standardised in preparation 
for cluster analysis to ensure equal prominence, as 
advised by Davis (1986). The clustering was carried 
out in two stages, in an approach similar to those 
of Fovell (1997) and Davis (1986) which allows 
atypical samples to be absorbed into local groups 
to maximise broadscale resolution. The clustering 
method chosen was the KMeans Divisive Clustering 
described by Hartigan (1975), and Hartigan & Wong 
(1979). This is an iterative partitioning procedure 
that assigns samples to a specified number of non- 
overlapping clusters. Although this procedure 
does not produce the usual dendrogram, it can 
accommodate large volumes of data and is resistant 
to outliers. Dissimilarity between the clusters was 
estimated using the Euclidean metric technique 
(Steinhorst & Williams 1985). 
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TABLE 1: Component loadings in %eq/L (percentage of the total equivalents per litre) as well as percent of variance explained 
for the first three Principal Components. 





The parameters used for the first stage of clustering 
are the percentages of the major cations, Na, Mg and 
Ca, and anions Cl, SO, HCO, and NO,, expressed as 
milliequivalent per litre (meq/L). These parameters are a 
common choice for classifying natural water chemistry, 
as illustrated by, for instance by Pacheco (1998), and 
Ground & Groeger (1994). For this procedure, the 
combined ground and surface water datasets, sorted 
by EC, were split into sets of 100 samples, each 
representing a narrow salinity range. Each set was then 
clustered into five groups, and these initial five groups 
were replaced by simulated samples based on the group 
medians. The new data set, consisting of 910 simulated 
samples, was then re-clustered into ten regional water 
types. Both the simulated data, and the original data 
that they represented, were labelled with the water type 
produced by the second stage of clustering. 


Because multi-stage clustering can produce distortions 
(Fovell 1997), the regional water types were further 
confirmed by a Principal Component Analysis (PCA) 
as described in McNeil, Cox & Preda (2005). As stated 
by anumber of authors, for instance Zuur et al. (2010), 
Jolliffe (2002) and Dunteman (1989), PCA makes no 
assumptions of normality, and is frequently combined 
with cluster analysis in water quality assessments (e.g. 
van Tonder & Hodgson 1986). It is useful as a visual 
indication of whether the clusters are sufficiently 
compact, whether they substantially overlap, or 
whether they contain natural separations. The input to 
the PCA was the 910 simulated samples produced by 
the first stage clustering and each labelled by the final 
ten water types produced by the second stage. Table 
1 lists the factor loadings on each of the first three 
principal components which contributed 94 % of the 
total variance. The table also includes the individual 
contributions of factors to the total variance. 


Figure 2 is a two dimensional representation of the 
factor loadings on the first two principal components, 
with the major ions for each final water type. This 
figure indicates that the separation, cohesion and 
compactness of each of the water types are sufficient 
to justify a separate identity. The strongest trend 
shown by the PCA is the inverse relationship between 
chloride and bicarbonate, with types I and F having 
high proportions of chloride. This is not surprising, 
since these are highly saline types, where enrichment 
in the conservative chloride ion could be expected. 
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FIGURE 2: Principal components analysis to verify water 
type groupings. The ten clusters are distinct from each other 
and sufficiently compact not to require further subdivision. 
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FGURE 3: Ionic distributions in terms of milliequivalents per litre for the water types comprising, showing the two main 


sequences, with seawater for comparison 


Of more interest is the inverse relationship between 
the sodium ion, and calcium, magnesium and sulphate 
collectively, with the more sodic types being F, G 
and D which cover a wide range of salinities. The 
third component differentiates the high proportion of 
sulphate in Type E water. 


For each water type, the variances of the major 
ions within the original dataset were plotted as box 
and whisker arrays similar to Schoeller diagrams as 
depicted in Hem (1989). These produced distinctive 
patterns which, together with the PCA outputs, enabled 
the water types to be grouped into broader sequences. 
Two major chemical sequences were defined, each 
covering a wide range of salinity. The Schoeller style 
diagrams for these sequences are shown on Figure 3, 
accompanied by that of typical seawater as presented 
by Hem (1989). Both sequences are widespread in 
occurrence, probably as a consequence of the size and 
complexity of the catchment. 


The first sequence, comprising water types A, C, B and 
I in increasing order of salinity, is titled the ‘alluvial 
sequence’ because it dominates the major alluvial areas 
as well as most surface waters. The cations in these 
water types are approximately equal in proportion (or 
balanced), with each of their medians representing at 
least 25% of the total. HCO, is the major anion, except 
in the highest salinity ranges, where NaCl begins to 
dominate. Previous studies, for instance McNeil, 
Cox & Preda (2005) and McNeil (1983), have shown 
this trend to be typical of both surface waters and 
alluvial aquifers associated with Palaeozoic deposits 
in Queensland. The second sequence, designated the 
“sodic sequence’, comprises water types D, G and F, 
is dominated throughout by Na and Cl ions, although 
the proportion of HCO, may be comparatively high at 
the lowest salinities. Figure 3 shows that this sequence 
evolves to resemble a marine ionic balance as salinity 
increases. The chemical characteristics of the water 
types are summarised on Table 2. 
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TABLE 2: Fitzroy groundwater chemistry types as derived from cluster analysis. Major ions are cations Na, Mg and Ca, and 
anions Cl, SO,, HCO, and NO,. 


Salinity is classed as very low if the median EC is below 200, low if between 200 and 500, moderate if up to 1500, high up to 
5000, and very high if in excess of 5000. It is classed as variable if the EC range is more than twice the median. 


A single dominant cation is named if its median exceeds each of other two by at least 25%, (or 30% in the case of sodium). 
Balanced major cations means that none constitutes less than 25% of the total in a median sample. 


Balanced Predominant surface water type in Fitzroy basin as a whole, 
and in most sub-basins other than those with extensive 
alluvium. Rare in groundwater apart from lower salinity sub- 
basins (Isaac, Comet and Nogoa). 
Balanced HCO,, Cl Low- Only common surface water type apart from A, and dominant 
Moderate in alluvial sub-basins (Lower Fitzroy, Isaac and Callide). 

Also most common shallow groundwater type apart from B. 
Significant proportions of Type C in surface and groundwater, 
as occurs in the main alluvial sub-basins, indicates substantial 

interaction. 


Na>Mg>Ca Cl Moderate Most common in shallow groundwater except in Comet sub- 
(SO, above basin, rare in surface water apart from Callide sub-basin. 


_ 


Saline Saline groundwater which is enriched in Cl with respect to Na 
which may have been replaced by divalent cations. Significant 
only in deep groundwaters of the eastern sub-basins, including 

the Callide and lower Fitzroy. 

HCO,>Cl Low- Most common deep groundwater type overall, and in southern 

Moderate and western sub-basins. Occasionally significant in shallow 

groundwater. In comparison to Type B, it is enriched in 
HCO, and Na at the expense of Cl and the divalent cations 
respectively, with little overall change in salinity. 
CI>HCO, Low- Minor type, similar to Type F with which it is associated, but 
Moderate | less saline and evolved towards marine-like composition. Also 
more common in shallow groundwater. 


Saline Mainly deep groundwater type, evolving towards a seawater 
composition. Most significant in Isaac and Mackenzie sub- 
basins. 
Cl, Net Very rare and insignificant low salinity type, high in sodium 
and nitrate, usually associated with irrigation areas. 


Balanced Moderate Rare, sulphate rich type, insignificant except occasional 
occurrences in shallow groundwaters of the Callide and 
Mackenzie, and historically in the surface waters of the Don 
and Dee Rivers. 


Moderate Rare magnesium bicarbonate dominated groundwater 
expected to be associated with basalts. 
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Variation of EC with depth 














Ec inuS/cm 
to 












































10 20 30 40 50 BO 70 
Depth inm 


Depth distribution of groundwater samples 


800 


























i*) 
Oo 
o 


400 


Number of samples 
Jeg jed uomodoig 








0 20 40 60 80. 
Depth in m 


FIGURE 4: Range of sampling depths for Fitzroy bores in the dataset, with median and 20-80 percentile range of EC to show 


variability with depth. 


In recognition of the vertical changes in chemical 
composition that commonly occur, both within 
and between aquifers, a boundary was required to 
differentiate shallow and deep groundwaters so as to 
confine zonal ranges within meaningful limits. For 
this broadscale study it was necessary to apply a single 
depth subdivision to the whole basin, maintaining the 
option to refine depth ranges for individual zones if 
necessary, during later investigations. Figure 4 shows 
that most of the samples were collected at depths of less 
than 30m, which is the maximum depth of most of the 
alluvial deposits. The depth distribution also appears 
skewed, with a tail of deeper bores which could be 
presumed to represent a range of aquifer types. In 
addition, the distribution of EC with depth, also shown 
on Figure 4, indicates a change in profile between 30 
and 40m. Above this level, the EC is declining and 
becoming less variable with depth, whereas at deeper 
levels its variability rapidly increases. On this basis, 
a depth of 30m was selected to regionally subdivide 
shallow and deep groundwaters. 


DEFINITION OF GROUNDWATER ZONES 

The regional groundwater zone boundaries were 
established on the basis of similar groundwater 
chemistry, with boundaries refined on geological and 
catchment divisions and in accordance with rainfall 
patterns. After the samples had been classified by 
water type and depth category, they were plotted as 


GIS layers (Figure 5) with bores displayed as points 
and surface water as pie diagrams (Figure 9) showing 
predominance of type. This enabled their geographical 
distributions to be superimposed on the environmental 
layers for visual comparison. 


The geological map layer (Figure 6) was drawn 
from the Queensland Department of Environment 
and Resource Management (DERM) SIR-QRY 
database, which in turn was extracted from the 
Rock Units Table held in the DERM MERLIN 
Database, founded on data from GSQ 1975. The 
geology of the Fitzroy basin has been described by 
a number of authors, for instance, Douglas et al. 
(2006), and Speck et al. (1968). The oldest rocks 
are metamorphics, sediments and volcanics, mostly 
of Palaeozoic age or earlier, which occupy the 
westernmost portion of the catchment and dominate 
the easternmost area which constitutes part of the 
Great Dividing Range. These deposits form hilly 
country, containing actively weathering minerals 
which influence groundwater either through direct 
contact or when redeposited as alluvium. These 
rocks were intruded late in the Mesozoic, and 
plutonic outcrops now form rugged landscapes with 
a local relief of over 300m, particularly along the 
eastern border of the catchment. The plutonic rocks 
may contain fracture zones, but are not normally 
significant in terms of groundwater. 
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FIGURE 5: Shallow and deep groundwater chemistry types in the Fitzroy. 
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FIGURE 6: Dominant geology in the Fitzroy. (DNRM ). 


State surface geology 1975 - Queensland. 2011-11-22 (publication). The polygons in this dataset are a digital representation 
of the distribution or extent of geological units within the area. Polygons have a range of attributes including unit name, age, 
lithological description and an abbreviated symbol for use in labelling the polygons. These have been extracted from the Rock 
Units Table held in the Department of Natural Resources and Mines MERLIN Database. The lines in this dataset are a digital 


representation of the position of the boundaries of geological units and other linear administrative features. The lines are 
attributed with a description of the type of line represented. 
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During the Mesozoic, and continuing into the Tertiary, 
thick layers of terrestrial sediment were laid down in 
the Bowen Basin, and also in the Surat Basin, which 
covers the southern catchment where it provides part 
of the southern recharge of the GAB, particularly 
in the headwaters of the Dawson. The sediments 
are mainly sandstone but also contain mudstone, 
volcanics, limestone, oil shale, coal and lignite. They 
form plains or dissected tablelands, with a local relief 
of 30 to 150m. Spring lines fed by underlying GAB 
aquifers in the south can contain unique ecosystems 
(Fensham et al. 2007). 


There were several phases of volcanism in the early to 
mid Tertiary, resulting in extensive basalt flows which 
remain in the west, particularly in the Comet and Nogoa 
subcatchments including Theresa Creek. Numerous 
basalt remnants are present in other parts of the Fitzroy, 
including the headwaters of both the Isaac River, and 
Callide Creek, upstream of Biloela in the Dawson 
subcatchment. Basalt contains important aquifers, and 
has been shown through other studies (McNeil 2002, 
McNeil, Cox & Preda (2005) to be associated with a 
water chemistry which is unusually high in magnesium. 
Later in the Tertiary, the landscape was affected by deep 
weathering, with the formation of duricrust and lateritic 
plains. Much of the present landscape is dominated by 
broad Quaternary flood plains where the weathered 
alluvium or colluvium has a local relief of less than 15m 
and is drained by deeply incised, low gradient streams. 
These older alluvia have been known to contaminate 
adjoining aquifers with saline seepages (McNeil 2002). 
Recent alluvium has built up within the older deposits 
around major streams, and these now constitute the 
major aquifers of the basin. 


More than 100 rock units were found to be present 
in the catchment, but the map on Figure 6 shows a 
simplified version of the surface geology in which these 
rock units have been divided into categories which 
have been shown in previous studies (for instance 
McNeil 1983 and McNeil; Cox & Preda 2005) to 
represent similar effects on bulk water chemistry. The 
geological categories were based on texture, mineral 
content, and ability to accumulate and release salt. As 
can be seen on Figure 5, there is significant alignment 
between the groundwater chemistry sequences. The 
sodic sequence tended to occur in the vicinity of 
Mesozoic sandstones, including material or recharge 
derived from them through erosion, transportation, 
and runoff. On the other hand, the alluvial sequence 
was generally aligned with alluvial deposits derived 


from the mixed sediments and volcanics of the Great 
Dividing Range, and/or from basaltic rocks. 


Long term average annual rainfall data sets were obtained 
from the SILO Australian climate archive, which stores 
daily time step rainfall interpolations in 5km grids (for 
interpolation methods, see Jeffrey et al. 2001). The data 
was processed for visualisation in ArcGIS 9.3.1 (Figure 
7). The effect of rainfall was aligned with other factors, 
such as topography and geology, and particularly with 
relative recharge rates available from the DERM SIR- 
QRY database, having been produced by departmental 
groundwater models. Figure 7 shows that high rainfall 
and recharge areas occurring near the coast, and along 
the ranges extending across the large southern Dawson 
sub-catchment, were dominated by alluvial sequence 
chemistry types, as would be expected from the geology 
in those areas. Sodic types, particularly those of higher 
salinity, tended to coincide with low rainfall and recharge 
areas. This led to some zones cutting across sub- 
catchment boundaries due to the observed continuation 
of water chemistry which was consistent with average 
rainfall trend. 


A general alignment was expected between 
groundwater chemistry and land use, because both 
are aresult of geology, soils, topography and climate. 
It was also anticipated that local anomalies in 
groundwater type within an otherwise homogeneous 
region could be associated with the impacts of 
specific land uses. The anomalous types could then 
be excluded from the calculation of baseline ranges. 
Land use was accordingly mapped into discrete 
categories (Figure 8) following the method of Cogle, 
Carroll & Sherman (2006). The category choices 
reflect the types of landscape intervention occurring, 
as opposed to land use descriptions based on outputs 
(Cogle, Carroll & Sherman 2006). The state-wide 
land use data set was obtained from the Queensland 
Land Use Mapping Project (QLUMP) which uses the 
Australian Land Use and Management classification 
(ALUM) Version 5 Feb 2002. The aggregated land 
use data set was then rasterised into 100m grids in 
ArcGIS 9.3.1. Although observed associations were 
not strong, broad alignment were noted between 
the alluvial sequence and undeveloped or forested 
areas, and between the sodic sequence and grazing 
lands; however there were many exceptions. Both 
chemistry sequences occurred in cropping areas, for 
example, and these areas often reflected complex 
groundwater chemistry with some indications of 
surface water and groundwater interaction. 
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After the zones had been delineated (Figure 9), 
percentiles of salinity and major ions were calculated 
for deep and shallow groundwater within each zone 
(Raymond & McNeil 2011). Shallow groundwater 
chemistry was also compared with that of the surface 
water within the same zone to identify possible areas 
of interaction. Some areas remained uncategorised, 
either because they contained no groundwater data, or 
because it was too scattered and variable to represent 
consistent groundwater systems. 


RESULTS AND DISCUSSION 

Despite the uneven distribution and density of 
groundwater data, 44 discrete chemistry zones were 
delineated across the Fitzroy Basin (Figure 9 and 
Table 3). Certain areas, such as parts of the northern 
and western fringes of the basin, have not been zoned 
because the geology and sparseness of data indicated 
that they do not contain substantial groundwater. 


Most classifiable areas of the basin were found to 
be dominated by the alluvial or sodic sequences. 
The alluvial sequence water types have no dominant 
cations, although the contribution of sodium increases 
with salinity, as does the proportion of chloride in 
comparison to bicarbonate. This sequence is generally 
aligned with Palaeozoic and basaltic rocks on the 
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FIGURE 7: Average annual rainfall (mm) from 1900 to 1999 
in the Fitzroy. (SILO). 


Total annual rainfall data sets and long term averages were 
obtained from the SILO Australian climate archive (data 
courtesy of Queensland Government http://www.longpaddock. 
qld.gov.au/silo/ and the Australian Bureau of Meteorology), 
which stores daily time step rainfall interpolations in 5km 
grids (Jeffrey et al., 2001) based on point rainfall data provided 
by the Bureau of Meteorology. 


eastern and western sides of the basin where rainfall and 
groundwater recharge are also relatively high. It is also 
found in the surface water and most alluvial deposits. 
As explained by McMahon, Bohlke & Lehman (2004) 
and Hem (1989), this chemistry is typical of actively 
weathering landscapes, which would include basalt 
terrains, as well as the mixed sediments and volcanics 
of the Great Dividing Range, and their associated 
alluvials. Such chemistry is consistent with base-cation 
release through chemical weathering of diverse mineral 
assemblages in subtropical settings with summer 
rainfall as described in Millot (1970), McNeil (1983), 
White and Blum (1995) and Murdoch et al. (2000). 


Conversely, the sodic water types align mostly with 
the granites, and with the terrestrial Tertiary and 
Mesozoic quartzose sandstones and their downstream 
alluvials, mainly in the lower rainfall areas of 
the central and southern parts of the basin. These 
sediments are related to the underlying Bowen and 
Surat Basins, which have relatively low recharge 
rates, and include some GAB springs. Authors such 
as Chae et al. (2006), Eriksson (1985), Kerala et al. 
(2007), Pinol & Avila (1992) and Shaji, Bindu & 
Thambi (2007) have attributed sodium enrichment in 
coarser rock types which are relatively low in ferro- 
magnesian minerals, to feldspar weathering, cation 
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FIGURE 8: Land use in the Fitzroy. (QLUMP). 
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FIGURE 9: Groundwater zones described in Table 3 with surface water types in the Fitzroy. The surface waters were collected 
as time series at gauging stations and are represented as pie diagrams indicating frequency of type. 
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TABLE 3: Summary of groundwater chemistry zones in the Fitzroy Basin. 


Zone 
- Dennison 


2 - Anakie- 
Sapphire 


3 - Mt Aldis- 
Concilliation 
4 - Blair Athol- 
Capella 


5 - Lower- 
Stephens 


6 - Krombit- 
Kariboe 


7 - Zamie- zane 


8 - Mt ef Soe 
9 - Sandy 


10 - Bee- 
Cooper 


11 - t-te | 


Isaac 


13 - Phillips- 
Fairbairn- 
Meteor 


14 - Prospect- 
Fitzroy 


Groundwater data 


Shallow Limited to 
moderate 
Deep Limited to 
abundant 
Shallow Limited to 
abundant 


Dominant ions 


Cations 


FICO, 


Balanced 


Balanced 


Balanced HCO 


Balanced HCO, 


HCO, 
HCO 
HCO 


3 


a 


Balanced 
Balanced 4 


Balanced 


Balanced 


EC (uScm-1) 
Percentile value 
80th 


Deep Limited 


Shallow Moderate to 


abundant 


zane Limited 


1,293 


Shallow st oe, — a Ca 


Suan Limited to 
moderate 


Shallow Moderate 


ps0] 1.200 | ,200 


Shallow Moderate to Sef co, HCO, oo Ea 
abundant 


(Na | CHICO, | HCO, 


Deep Moderate to 


abundant 


761 832 
fer — 


Shallow | Moderate to O, 1,150 2,509 
abundant 

Deep Moderate to Na HCO, 756 1,420 2,150 
abundant 


Shallow CI,HCO, | 1,006 1,619 2.765 
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Zone 


15 - Upper- 
Nogoa 


16 - Thangool- 


Jambin 
17 - Lotus- 
Connors- 


Funnel 


18 - Rannes 


19 - Dee-Don 


20 - Robinson- 


Baffle 


23 - Dysart 


24 - 
Tomahawk 


25 - Lower- 
Meteor 


26 - Wharton- 
Claude 
27 - Cockatoo 


USING STATISTICAL AND CONCEPTUAL METHODS 


Depth Groundwater data | Dominant ions EC (uScm-1) 
Percentile value 
Cations | 20th 50th 80th 
Deep Limited to 1,340 
moderate 
Shallow | Limited }229 [1,050 050 Masts 515 
7 —_} 1.050 {ee 


Deep Moderate to 2,721 3,900 7,200 
Fie laine 


moderate 


} Shallow |Limited | 
Soon ee to uno, THEO: 2 2,060 
moderate 
o 


Shallow | Limited to Cl, HCO, 1,033 
moderate 


Deep Limited to 
abundant 

Shallow | Moderate to Balanced Cl 1403 2,220 3,722 
abundant 

Deep Limited to 2,496 3,465 7,450 
moderate 


Deep 
Shallow 


Limited to 
moderate 


Abundant 
Limited 


Deep 


Shallow 


Moderate 


Limited to 
moderate 


Moderate to 
abundant 


Deep 


Shallow Moderate to 


abundant 





Deep Moderate 


Shallow | Limited 
Deep Abundant 


Shallow Moderate 





54 PROCEEDINGS OF THE ROYAL SOCIETY OF QUEENSLAND 


Zone Depth Groundwater data | Dominant ions EC (uScm-1) 
Percentile value 
Deep Limited to Na 
28 - fen dy abundant 
Conciliation- = 
Mimosa Shallow | Limited to Na HCO, 300 615 1,644 
abundant 
890 
530 
536 
2,151 


29 -Injune 689 1200 | 2.205 
Palmiree | Shallow | Moderate [Na | CLHCO,_| 


Deep Limited to Na HCO, 
30 - Retro- moderate 
CO, 


2,515 


Darrylar Shallow | Limited to Na H 
moderate 


Deep Limited to Cl, HCO, 3,150 3,540 
moderate 
31 - Sandhurst 


2,050 2,050 | 5 165 
}1440 | 440 6,040 


6,3 aa 8,712 
Castle Shallow | Limited to Na 3,700 2,158 
moderate 
Deep Limited to 3,419 6,100 16,000 
34 - Isaac- moderate 
Dawson Shallow | Limited to Na Cl 498 2,150 8910 
moderate 
Deep Limited to 4,103 13,604 
35 - Brown- aaa seit 


Humbolt 
Blackwater ee st fs fe Cl, HCO, ee 100 ae — 


| Deep 


37 - Clermont- 
Theresa Shallow en to a 400 441 2 
moderate 
}Deep |Deep | Limited | [Na | CLHCO, | HCO, 1333 | 333 1675 ,675 }1675f2.450 | 
38 - Callide Shallow pie _ to ——— 1,050 1,610 
moderate 


|Deep —_| |Deep | Limited }1270 [1,422 | 422 


39 - 
Marlborough Shallow = to eS _= 314 ae 646 ie 725 
moderate 


|Deep | |Deep | Limited | Balanced 805 1363 ,363 1616 | ,616 


40 - Balm 
*_[ Shallow ese ee oe 


Shallow Limited to Na Cl, HCO, 2.450 4.200 
moderate 

Deep Limited to Na 
33 - Zamie- moderate 
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Zone Depth Groundwater data | Dominant ions EC (uScm-1) 
Percentile value 


rer 


42 - Black 
Hole 


Comet 


44 - Dee-Don- 
Callide 


exchange with clays, and a reduction in the proportion 
of calcium in relation to other cations through calcite 
precipitation in a drier environment. Some water types 
occur which are outside the two main sequences, but 
these are rare and related to local factors. 


Significant differences in bulk chemistry between 
deep and shallow groundwater may indicate overlying 
but unconnected groundwater systems. About half of 
the 44 identified zones reflected similar bulk chemistry 
in both deep and shallow groundwater, about a 
third reflected significant differences in salinity and 
chemistry across depths, and the remainder were 
lacking either deep or shallow groundwater data to 
be assessed. Differences in chemistry across depths 
were mostly found in alluvial sequence zones, 
where cations were more balanced in the shallow 
groundwater. However in most of these zones, apart 
from Zones 2 and 27, there was limited data for either 
the shallow or deep groundwater, or both, reducing 
overall confidence in these observations. 


Groundwater and surface water systems are 
considered unlikely to be physically connected where 
their proportions of major ions are substantially 
different. On the other hand, there is a greater 
likelihood of surface water and groundwater being 
in recent contact, either through recharge or baseflow 
if their bulk chemistries are reasonably similar. The 
distribution of chemistry types between surface water 
and shallow groundwater was examined in each zone 
to determine the likelihood of significant interchange 


Shallow Limited to N 
moderate 


Deep Limited to Balanced 
moderate 

Shallow | Limited to Balanced C 1 
moderate 


a 
Deep Limited to Balanced 
43 - Lower- moderate 
| : et 





z 


[500 [a 


(Figure 9). The alluvial chemistry sequence usually 
extended into the adjacent surface water, albeit with 
a reduced salinity. The surface water then tended to 
maintain its character downstream, even extending 
over sodic groundwater sequences which are very 
rarely expressed in the surface water. The surface 
waters were found to be mostly fresh and belonging 
to the alluvial sequence (mostly water type A). 
Potentially significant interaction was noted in some 
areas where groundwater chemistry was reflected at 
the surface, particularly around Mt Morgan (Zone 41). 


8) 
Cl 
Cl 
Cl 
Cl 

l 


In terms of reliability, it is conceded that the results 
of a broad scale desktop study such as this, based on 
historic data, has limitations in terms of spatial and 
vertical precision. There area variety of aquifer systems 
which may or may not interact, and broad areas where 
groundwater is only found sporadically. Ground and 
surface water sampling although extensive, has also 
been largely unsystematic, with records collected 
from non-uniform locations in the landscape, and 
with variable quality control. The spatial distribution 
of data is therefore uneven, with most emphasis 
having been given to irrigation areas. Subsequently, 
the zones were delineated and characterised with a 
high degree of confidence where groundwater data 
were abundant and relatively homogenous in water 
type. However, less certainty was possible in areas 
where the data were inadequate or biased, or where 
the water chemistry patterns were very complex. Such 
complexity is common in areas with small isolated 
bores, with limited storages, that are likely to differ 
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widely in chemistry because of local environmental 
and land use influences. In particular, it would be 
desirable to carry out more detailed investigations 
in Zones 2, 3, 13, and 35 as the data indicate strong 
vertical zonation which may not align with the 30m 
regional boundary, as well as high spatial variability. 
Despite these reliability concerns, all defined zones 
reflect a relatively consistent chemistry, and are 
sufficiently aligned with relevant geological, climatic 
and land use indicators to enable an association with a 
particular water type for a given area. 


In recognition of the heterogeneous nature of the 
data, and the fact that sampling procedures were 
not documented during the period when most of 
them were collected and analysed, validated tests 
were carried out to identify unreliable samples, and 
robust, nonparametric multivariate procedures (fully 
described in NcNeil 2005) were used for statistical 
assessment. Shand et al. (2007) and Edmonds et al. 
(2003) emphasise that characterising unsystematic 
water chemistry data through a purely statistical 
approach can produce erroneous results, and that 
conceptual information should be incorporated in 
the assessments. This approach requires subjective 
judgement in delineating zonal boundaries, and in 
defining the degree of subdivision and levels of 
uncertainty varies across the basin. In terms of the 
confidence associated with parameter ranges, Shand 
et al. (2007) state that range tables for baseline water 
quality should ideally represent conditions that 
preceded all development. However, data for these 
times are rarely available, as is the case with this 
study. The estimated baseline ranges can therefore 
only reflect the existing condition of the groundwater 
as discernible from the historic data, and cannot 
account for changes which may have occurred in the 
time before it was collected. Baselines ranges which 
reliably reflect those of pre-development ranges 
could only be achieved through resource-intensive 
modelling, focussing on areas where change is likely 
to have been significant, however this is beyond the 
scope of this study. 


CONCLUSION 
This study presents a methodology to define reasonably 
homogeneous zones of groundwater chemistry, 
when only historical data is available, and this may 
be unevenly distributed and confined to salinity and 
major ions. Such a methodology is needed to produce 
baseline parameter values to fulfil the requirements 
of the precautionary principle to guide adaptive 


management pending more intensive investigations. 
The methodology is based on multivariate analysis 
of water chemistry data, combined with conceptual 
information supplied as multiple GIS layers. It was 
applied to the Fitzroy Basin, a large catchment with 
complex groundwater. Baseline water chemistry 
ranges are developed for both shallow and deep 
groundwaters within each zone. The addition of 
surface water samples to the groundwater data for 
combined multivariate analysis is introduced to 
identify areas of interaction. 


The nature and distribution of groundwater chemistry 
in the basin has been mapped, and baseline water 
quality ranges produced for reference. In particular, 
broad alignments have been defined between two 
sequences of water type, each evolving through a wide 
range of salinities. The zonation of the sequences with 
their respective geological and climatic indicators 
held true for most areas, and was useful in justifying 
the classification of dominant groundwater chemistry 
in areas of less certainty. 


Continuous monitoring is considered an essential 
element of groundwater quality assessment programs 
by many authors, for instance, Shand et al. (2007) 
and Nieto et al. (2005). Further monitoring and data 
recording is therefore recommended to clarify the 
dominant groundwater chemistry in complex areas or 
those with sparse data. This would clarify the risks to 
local ecosystems, and enable the detection of long term 
trends in areas where groundwater and surface water 
interaction may be associated with human activities. 
There is also the potential to further develop the 
methodology through more computational approaches 
to GIS mapping. 
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WHY THERE IS NO GENETICS IN PHYLOGENETICS: MISCONCEPTIONS, 
CLARIFICATION AND SOLUTIONS 


SCHARASCHKIN, T 


Fundamental misconceptions regarding some basic phylogenetic terminology are presented in this 
opinion piece. An attempt is made to point out why these misconceptions exist and what may be caus- 
ing the misapplication of terminology. Clarification is providing via basic definitions and simple ex- 
planations. Differences between the scientific fields of genetics and population genetics are discussed. 
The appropriate use of terminology is advocated and alternative terms are proposed to eliminate one 
potential source of confusion. It is suggested we use ‘sequence data’ instead of molecular data and 
“non-sequence data’ instead of morphological data in the field of phylogenetics and systematics. 


Scharaschkin, T., Earth, Environment and Biological Sciences, Science and Engineering Faculty, 
Queensland University of Technology, Brisbane, Queensland, 4001, Australia. 


INTRODUCTION 

As a plant systematist, actively involved in 
research and teaching, I interact with a wide range 
of students and professionals in the area of plant 
sciences. I have frequently encountered fundamental 
misconceptions regarding the meaning of the terms 
Cladistics, phylogenetics, taxonomy and systematics. 
Subsequent consequences of misconceptions include 
a lack of understanding for the research methods 
involved and the interpretation and implications of the 
results. These misconceptions are not only prevalent 
in the undergraduate and postgraduate community, 
but also surprisingly encountered in conversations 
with academics and researchers. This opinion piece is 
an attempt to point out some of the misconceptions 
whilst offering an explanation as to why they exist and 
how we can avoid misapplication of these terms. 


COMMON MISCONCEPTIONS 

CLADISTICS AND PHYLOGENETICS 

The terms cladistics and phylogenetics are often 
considered (incorrectly) to differ in the type of data 
used: morphological data (cladistics) and genetic 
data (phylogenetics). This leads some people to think 
that cladistics is old fashioned while phylogenetics 
is a more modern technique. Students often express 
the view that phylogenetics is a better approach 
than cladistics as problems of assessing homology 
can be avoided. On the other hand, more seasoned 
researchers have sometimes suggested that cladistics 
requires a thorough understanding of the biology, 
especially morphology, of the research organism, 
while phylogenetics is the quick and easy lab approach 
that can be conducted without much knowledge of the 
biology of the organism of interest. 


TAXONOMY AND SYSTEMATICS 

The overall use of these two terms is less frequently 
encountered, possibly as a consequence of less 
time being devoted to these topics in undergraduate 
teaching or actively pursued in research programs. 
Undergraduate and postgraduate students often 
acknowledge having heard these terms, but if probed 
further, the lack of clarity becomes apparent. When 
asked what the similarities or differences are, the 
hesitant but most popular answer provided is that 
taxonomy utilises morphological data and/or is old- 
fashioned while systematics relies on genetic data 
and/or is a modern method. 


MORPHOLOGICAL AND MOLECULAR DATA 
The type of data used in phylogenetic research is 
frequently misunderstood. The most frequently 
encountered misconception 1s that morphological data 
is restricted to the morphology of the organism while 
molecular data involves genes. And as mentioned 
above, there is also the accompanying perception that 
cladistics and taxonomy utilise morphological data 
while phylogenetics uses only molecular data. 


CLARIFICATION 
Many of the misconceptions outlined above probably 
stem from a lack of formal introduction to the 
specific subject. Without going into the details of the 
development of the field of phylogenetics and the 
history of key researchers who have influenced this 
science, some of these misconceptions can be easily 
cleared up by simply paying attention to etymology and 
definitions. I have relied on a number of resources (cited 
only once and collectively) to present the information 
in a greatly simplified form (Felsenstein, 2004; Hillis et 
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al., 1999; Hillis & Wiens, 2000; Kitching et al., 1998; 
Wiens, 2000; Wikipedia). Definitions are summarised 
in Table | to serve as a quick reference guide for future 
use by readers. In a few cases more than one definition 
has been provided to facilitate a more complete or 
formal understanding of the terminology. 


CLADISTICS AND PHYLOGENETICS 

The origin of the word clade comes from the Greek 
word k/ados (meaning a shoot or branch). The 
term clade refers to “a group or subset” and can be 
applied to any collection of objects, be they animate 
or inanimate. The term cladistics simply refers to the 
methods, including criteria, used to determine clades. 
Clades can be represented diagrammatically, such 
as in the form of a branching tree. The diagram of 
relationships is called a cladogram. In evolutionary 
biology, the criterion applied for grouping organisms 
together is the presence of shared evolutionary history 
or descent from a common ancestor. 


The term phylogenetics is derived from the Greek 
terms phylon (meaning tribe or clan) and genetikos, 
the adjective of the word genesis (meaning origin or 
birth). Phylogeny therefore refers to the relationships 
among groups. The term phylogenetics refers to the 
methods used to reconstruct a phylogeny. The diagram 
of relationships is called a phylogenetic tree. In 
biology, a phylogeny is the evolutionary history of a 
group. As this is something that happened in the past, 
it is unlikely we will ever know the true phylogeny of 
the group we are interested in, therefore the best we 
can do is to infer the phylogeny. As with cladistics, 
the criterion applied for grouping organisms together 
is the presence of a shared evolutionary history L.e., 
descent from a common ancestor. The hypothesis of the 
evolutionary relationships of a group can be represented 
diagrammatically, in the form of a branching tree. 


It is worth noting that this explanation of the words 
cladistics and phylogenetics does not require any 
information on the type of data or analytical method. In 
biology, the two terms are therefore synonyms and can 
be used interchangeably as long as the basis for grouping 
organisms assumes a shared evolutionary history. 


TAXONOMY AND SYSTEMATICS 

The term taxonomy encompasses the field of science 
that deals with the description, identification, 
nomenclature and classification of organisms. There 
are no firm rules about how to classify or group 
organisms. Historically, much of the classification has 


relied on the opinion of experts (intuitive taxonomy) or 
overall similarity (phenetics). Systematics is the field 
of science that deals with the description, identification, 
nomenclature and classification of organisms in an 
evolutionary context. In other words, it incorporates all 
aspects of taxonomy, but uses the shared evolutionary 
history of organisms as the criterion for classification. 


As with the distinction between cladistics and 
phylogenetics, the difference between taxonomy and 
systematics does not depend on the type of data used 
but on the philosophical question of how to group 
organisms so as to provide as reliable classification. 


MORPHOLOGICAL AND MOLECULAR DATA 
The cause for this last misconception is probably more 
the fault of the practitioners involved rather than the 
observers. In practice, the term ‘morphological data’ 
encompasses a wide range of information, including 
data obtained by studying the morphology of 
organisms but also anatomical, chemical, behavioural, 
ecological and other forms of data. On the other hand, 
the term ‘molecular data’ consists of DNA sequences, 
which may or may not code for actual genes. For those 
seeking a detailed explanation of the use of molecular 
or morphological data and the inherent problems and 
underlying assumptions of either type of data, the 
book “Phylogenetic Analysis of Morphological Data” 
(Wiens, 2000) is recommended. The title itself should 
be a cause for concern if the reader still believes that 
phylogenetics only applies to the use of genetic data! 


A common theme regarding the misconceptions 
presented here seems to lie with the misunderstanding 
that the terms taxonomy, systematics, cladistics and 
phylogenetics depend on the type of data used by 
evolutionary biologists. An additional source of 
confusion could be the fact that the term genetics 
shares part of the same root as that for phylogenetics, 
i.e. the part to do with genesis (= origin or birth). In 
biological sciences, genetics is the field that deals 
with the transmission and variation of inherited 
characters (i.e., chromosomes and DNA). It does 
not seek to reconstruct evolutionary relationships. A 
third potential source of confusion is the incorrect 
application of the term genetics to the field called 
population genetics and the related field of population 
biology. Population genetics is the study, at the 
population level, of allele frequencies and change 
under the influence of evolutionary processes such as 
natural selection, genetic drift etc. Population biology 
(often used interchangeably with population ecology) 
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Table 1. Definitions of terms presented in the text. 


Term 
Clade 
Cladistics 


Cladogram 


Genetics 
chromosomes and DNA). 


Phylogeny 
Phylogeny 
Phylogenetics 
Phylogenetic tree 
Population biology 


Population genetics 


A group or subset 


Methods used to determine clades 
Methods of classification that groups taxa hierarchically into discrete sets and subsets 
A diagram of relationships 


The field of science that deals with the transmission and variation of inherited characters (.e., 


The historical relationships among lineages of organisms or their parts (e.g., genes) 
The evolutionary history of a group 


Methods used to reconstruct the phylogeny 
A diagram of relationships 


The study of populations of organisms, especially the regulation of population size and various 
life history traits 


The study, at the population level, of allele frequencies and change under the influence of 


evolutionary processes such as natural selection, genetic drift etc. 


Systematics 


Taxonomy conducted in a phylogenetic framework 


The field of science encompassing description, identification, nomenclature and classification 
using shared evolutionary history (phylogeny) as the basis for the classification proposed 


Taxonomy 


is a Study of populations of organisms, especially the 
regulation of population size and various life history 
traits. Population genetics is much more closely allied 
with phylogenetics and systematics. The difference 
between population genetics and phylogenetics lies 
in the scale at which research questions are directed. 
Population genetics, as the name implies, is directed 
at the level of populations, where as phylogenetics 
(=cladistics) involves research at the species or higher 
taxonomic levels. 


POTENTIAL SOLUTIONS 

It is hoped that this overview, whilst based solely on 
anecdotal evidence of misconceptions, will highlight 
the need to be mindful of the proper use of terminology. 
If, as scientists, we do not use or apply terminology 
accurately, then our ability as communicators of 
science, especially to the broader public, is greatly 
damaged. We could be propagating misconceptions 
by simply truncating terminology. In using the term 
genetics, when we mean population genetics, we are 
conveying incorrect information, not only about the 
type of data but the entire field of research. 


If the terminology we are using is creating confusion, 
then we should provide clarification. But it also 
worth considering whether the problem lies with the 





The field of science encompassing description, identification, nomenclature and classification 


words we use, which might have a historical basis but 
have now outgrown their original usage. If the terms 
themselves are the cause of the problem, then it is 
up to the main users of that terminology to provide 
the solution. In the early days of phylogenetics, the 
only available source of data was morphology, but 
this is no longer the case. As a user of the terms that 
are the source of some confusion, I therefore propose 
the use of alternative words for the terms molecular 
and morphological data. I suggest we use “sequence 
data’ instead of molecular data and ‘non-sequence 
data’ instead of morphological data in the field of 
phylogenetics and systematics. 


In conclusion, by paying attention to the definitions, 
using formal and precise terminology and being aware 
of etymology of words, it becomes quite apparent that 
there is no genetics in phylogenetics. 
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